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abstract:   This  study  is  paralleled  after  J.  Hubbard's 
treatment  of  collective  modes  in  the  uniform  electron  gas 
and  differs  mainly  in  the  mathematical  complexity.   Large 
contributions  to  the  correlation  energy  are  identified  with 
collective  modes  and  this  identification  provides  conditions 
for  the  frequency  and  lifetime  of  these  modes.   The  contri- 
butions to  the  correlation  energy  are  analyzed  via  diagramatic 
techniques  where  the  Coulomb  interaction  between  the  electrons 
is  replaced  by  a  modified  interaction.   Only  first  order  or 
"bubble"  diagrams  are  considered  and  this  is  equivalent  to 
the  random  phase  approximation.   The  modified  interaction 
is  obtained  from  a  certain  integral  equation  whose  solution 
is  approximated  by  a  variational  principle.   The  kernel  of 
this  integral  equation  depends  upon  the  density  operator  for 
the  system  of  electrons  with  an  external  potential.   The 
solution  of  the  problem  is  formulated  in  two  different  ways, 
by  using  two  different  approximations  for  a  closed  form  of 
the  density  operator.   The  first  method  uses  a  variational 
principle  for  the  Green's  function  of  the  system  and  the 
second  method  uses  a  corrected  Thomas-Fermi  density  operator 
which  is  valid  for  slowly  varying  potentials.   Both  methods 
require  numerical  techniques  for  their  completion  and  these 
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are  carried  out  only  for  the  second  .nethod.   An  application 

is  made  to  atomic  type  systems  and  two  undamped  collective 

,  o 
modes  are  fo-ond.   One  is  a  very  singular  resonance  at  54A 

o       ^'^ 
and  the  other  is  a  broad   resonance  between  .24a  and  .46A. 
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CHAPTER  I 
Introduction 

It  is  known  that  a  free  electron  gas  exhibits  an 
undamped  collective  oscillation  at  the  well  established 
plasma  frequency.   We  would  like  to  study  how  an  external 
potential  affects  the  collective  oscillation.   In  order 
for  a  collective  mode  to  be  well  enough  defined  to  be 
observed  experimentally,  it  must  not  be  heavily  damped. 
This  Implies  that,  in  particular,  it  is  important  to  know 
how  the  inhomogenelty  of  the  system  affects  the  lifetime 
of  the  collective  mode. 

Several  theoretical  methods  of  studying  the  collective 
motions  of  a  many  body  system  exist.   For  example,  in  the 
theory  of  Tomonaga   a  canonical  transformation  is  made  so 
that  some  of  the  new  coordinates  are  directly  related  to  the 
collective  modes  of  motion  and  the  remaining  new  coordinates 

are  related  to  quasi-particle  modes.   In  the  theory  of 

2 
Bohm  and  Pines,   extra  coordinates  are  introduced  together 

with  an  equal  number  of  subsidiary  equations  so  that  the 

correct  number  of  degrees  of  freedom  is  preserved.   A 

canonical  transformation  is  then  made  so  that  the  extra 


Tomonaga,  S.,  Progr.  Theor.  Phys . ,  Osaka  1_3  (1955)  ^o?. 

2 

Bohm,  D.  and  Pines,  D.,  Phys.  Rev.  92  (1955)  509. 
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coordinates  when  transformed  are  related  to  the  collective 
modes  and  the  original  coordinates  whan   transformed  are 
associated  with  quasi-particle  modes.   In  "both  of  these 
methods  J  the  transformed  Hamiltonian  is  separable  with 
respect  to  the  collective  motion  part.   Though  both  of 
these  methods  are  quite  successful  when  studying  the  uniform 
electron  gas,  they  are  difficult  to  use  for  a  system  with 
an  arbitrary  external  potential.   First,  it  is  not  easy  to 
generate  the  above  mentioned  canonical  transformations  for 
an  arbitrary  system.   Second,  in  the  Bohm  and  Pines  theory 
it  is  difficult  to  find  eigenfunctions  which  satisfy  the 
subsidiary  equations.   Furthermore,  with  these  methods, 
the  interaction  between  the  collective  modes  and  the 
quasi-particle  modes  is  not  treated.   This  makes  the  damping 
problem  difficult  to  study. 

In  a  series  of  papers  in  1957    Hubbard  discusses 
a  method  which  is  applicable  to  an  arbitrary  fermion  system 
in  an  external  potential.   His  theory  is  based  upon  an 
infinite  perturbation  series  which  in  principle  is  an  exact 
solutioQ  of  the  many  body  problem  and  therefore  contains  all 
the  physical  effects  including  the  collective  motion. 
In  his  theory,  the  perturbation  series  is  simplified  by 


^  J.  Hubbard,  Proc.  Royal  Soc.  (London)  A24q  (1957)  539-560. 
J.  Hubbard,  Proc.  Royal  Soc.  (London)  A244  (1958)  199-211. 
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replacing  the  interaction  between  the  particles  by  a 
modified  interaction.   This  modified  interaction  is 
equivalent  to  the  modification  of  the  ordinary  interaction 
between  the  particles  by  the  polarization  of  the  medium 
represented  by  the  remaining  particles.   Thus^  Hubbard's 
theory  is  equivalent  to  a  dielectric  theory  viewpoint 
where  the  electron  gas  is  regarded  as  the  dielectric 
medium.   Using  Hubbard's  approach,  the  damping  problem 
is  easily  treated.   The  specific  calculations  done  in  this 
paper  are  based  upon  these  methods  of  Hubbard.   Therefore, 
in  the  next  section  we  present  a  detailed  su.'nmary  of 
Hubbard's  many  body  perturbation  theory. 
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CHAPTER  II 
Hub  1)3 rd  's  Many  Body  Perturbation  Theory-^ 
1 .   Introduation 

The  problem  is  to  determine  the  energy  spectrum  and 
wave  functions  of  a  gas  of  Permi-Dirac  particles  which 
interact  with  one  another  through  an  instantaneous  two-body 
potential  v.   It  is  assumed  that  the  particles  are  moving 
so  slowly  that  relativistic  effects  can  be  neglected. 
The  Hamiltonian  for  such  a  system  is 

H  =  Hq  +  h' 


where 


2 

Ho  -    2Z  a  +  1=  "(5?. )  +  YZ  V(x.  ) 

IX  1 


i7j 


h'  =  ZZ^(^i  -  ^i)  -  ^v(^.) 


2 

Pi 
5   —  is  the  kinetic  energy  of  the  particles. 

j7~  2m 


U^x.  )   is  the  potential  energy  of  each  particle  in 
an  external  field. 

^   v(x.  -  X.)  is  the  interaction  energy  of  the  particles 
and  V"(x.)  is  introduced  as  a  self-consistent  field.   The 


5 


J.  Hubbard,  0£.  Git. 
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need  for  V(x.)  as  well  as  a  defining  equation  for 
for  V(x..  )  will  be  discussed  in  more  detail  in 
Section  J   of  this  chapter. 
H'  is  treated  as  a  perturbation. 

2 .   Th_e_  Perturbation  Series. 

The  perturbation  H'  is  assumed  to  be  switched  on 
adiabatically  between  t  ^  -  oo  and  t  ^  0  and  adiabatically 
switched  off  between  t  =  0  and  t  =  oo .   Then  a  system  in  an 
eigenstate  ^„  of  H^^  at  t  =  -  oo  will  slowly  change  into  an 
eigenstate  Tp   of  H.   This  result  can  be  stated  in  the 
followina:  formi 


f   =   lim    3  (0,-aD)^Q  /  (^qIs  (0,-cx3)1^q) 
a  ->  +0 

is  an  eigenstate  of  H  with  the  energy 

E  =  Eq  +  AE  =  Eq  +  lim   (^qIh'  S  (O, -oo  )  1  ^q) /(^q  |  S  (O, -oo)|  V'q 

a  ^  +0 

where  E„  is  the  energy  of  the  anperturbed  state  f^ 
and  3  (t^t')  is  the  solution  of  the  equation 

(2.1)   ih  ^^-  S^(t,t')   =   H^(t)  S^(t,t') 

with  the  boundary  condition  S  (t',t')  =  1-,  where 

it  A  H^   ,   -it  A  H^   -a|tl 

^  T_r    /-s  ^  -, 


(2.2)   H  (t)  =  e       ^   H   e 
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The  limit  a  -^  +0  means  a  -♦  0  infinitely  slowly. 
Equation  (2.1)  is  equivalent  to  an  integral  equation 


(2.3) 


S  (t) 


1  + 


1 
ill  t 


H  (f  )  S  (f  )  dt' 


■CO 


where  S  (t)  stands  for  S  (t.-oo). 

The  perturbation  series  is  obtained  by  solving  (2.2) 

by  iteration. 


_oo 

n=^l 


t,    t   , 

1    n-1 


(2.4)   S^^(t)  =  1  +  ^   {^fl   dt^  J  dt^.../  dt^  -\(t3_)H^(t^ 


-oo    -oo    -oo 


...  H  (t  ) 
a^  n' 


oo 


n=l 


dt 


1 


dto  ... r  dt^  P[H  (t. )...H  (t^)] 
2    J    n   '■a^l'    a^n'-' 


00    -oo 


-co 


(P  is  the  chronological  ordering  operator  which  orders 
earliest  times  to  the  left.) 

Using  the  notation  of  field  theory,  the  perturbation 

H'  is 


2.5)   H 


1 
2  J 


[^(x'  )^(x'  )v(x  -X'  )^(x)^'(x)]  dx  dx' 


-  -^  N  v(0)  -  f  V'(x)  V(x)  V'(x)  dx 


where  v(x-x')  is  the  mutual  potential  energy  of  two  particles 
at  X  and  x'  and  N  i£  the  number  of  particles  in  the  system. 
The  second  term  subtracts  the  self-energies  of  the  particles 
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due  to  the  interaction  v  since  these. are  included  in  the 
first  term.   The  last  term  represents  the  contribution  from 
the  self-consistent  field.   ^(x)  and  ^(x)  are  the  particle 
field  operators  and  can  be  written  in  terms  of  an  ortho- 
normal  basis  [u.)  as 


1 


i 


->\  — 

i 


where  n.  and  n.  are  creation  and  destruction  operators  for 
the  particle  in  state  i  and  since  we  have  a  fermion  system 
n.  and  n.  obey  the   antieommutation  rules 


0:=[n^,n.]^  -    [n..nj]^,   t^^^jV=   ^i 


Here  the  u. (x)  will  be  taken  as  eigenfunctions  of  the 
unperturbed  one  particle  Hamiltonian, 

(2.6)   [|^  +  U(3?)  +  V(^)]u.(5?)  .  E.u.(5?)  . 

The  second  term  of  H  ,    i.e.  -  p-  N  v(0)^  is  a  constant 
and  will  be  included  from  now  on  in  H^  and  not  in  H  . 
If  we  write 

X   for   (xjt) 

V  (x-x')   for   v(x  -X' )5(t-t')e"^'*l        , 

o'  -i  —  E 

^(x)  =  ^(x,t)  =  Yll  ^iM'^i   =  IZ  Ui(^)e   ^   ""  n. 

i  i 

^(^)  =  y~  u^(x)n^ 


V(x)  =  V(x)e-°^l^ 


-7- 


and  substitute  equation  (2.5)  for  H   into  equation  (2.2) 
for  H  (t),  then  equation  (2.4)  for  S  (t)  tiecomes 


oo     -j^        -j^     n 


(2.7)  s^(t)  =^1  +|i:ht  (if) 


X  P{[|- J  dx^ax[{y^{x[)f{x[)v{x^-x[)^H^j^)-^{x^) 
-  J  dx^   ^(x-l)V(x^)^(x^)]    X   ... 


^    [|J  ^V<(^(<)^^<)^(^n-<)"^(^n)^(^n) 


-  Jdx^   ^(x^)V(x^)^(x^)]} 


where  the  integration  is  over  all  of  space  and  for  time 
from  t  =  -ootot=  t. 

5 •   Diagrammatic  Interpretation  of  the  Perturbation  Series 

It  is  convenient  to  re-express  the  particle  field 
operators  ^(x)  and  ^(x)  as 

^(x)  =  ^+(x)  +  ^~(x),  fix)   =   ^"^"(x)  +  ^"(x) 

where 

\inoc_c 
^   (^)    =        ^       u. (x)n.         (destroys   particles) 


1 

^  (x)  =   ^    u. (x)n.    (creates  holes) 

i 
occ 
V'  (^)  =   ^    u.(x)n.    (destroys  holes) 

i 
unoGC 

V  (x)  =   5ZI  ^- (^)^i    (creates  particles) 
i 
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We  are  interested  only  in  the  operation  of  S   on  a 
particular  eigenstate  f^   of  E^.      It  is  assumed  that  Tp^ 
is  a  non-degenerate  eigenstate  of  H„  in  which  certain 
definite  states  i  are  occupied.   The  state  ^„  is  then 
taken  as  the  "redefined  vacuum  state". 

We  then  have 

(2.8)  ^"(x)V^Q  =   0 

^~(x)^Q  -   0  . 

The  P  product  in  equation  (2.7)  is  resolved  in  such 
a  '/vay  as  to  move  all  the  operators  ^  ^  ^  to  the  left  and 
all  the  operators  f' ,    ip'    to  the  right,  making  use,  of 
course^  of  the  correct  corajmitation  rules.   After  this  has 
been  done,  because  of  equations  (2.8),  all  the  terms 
involving  Tp~   and  V'  can  be  omitted. 

It  is  then  found  that  the  remaining  terms  contributing 
to  S  (t)^„  can  be  conveniently  classified  according  to 
diagrams  similar  to  Feynmann  diagrams.   There  is  a  one  to 
one  correspondence  between  each  diagram  and  each  term 
contributing  to  S  (t).   We  nox^;  give  the  prescriptions  for 
drawing  these  diagrams  and  then  for  calculating  the 
corresponding  contribution  to  S  (t). 
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1.    Mark  p  points  on  the  diagram  and  label  them 

x-,jXp,...,x  .   Mark  another  set  of  p  points  and 
label  them  with  x-,  ^Xp,  .  .  .  _,x  .   Join  the  pairs  of 
points  x.;,x.  by  interaction  lines  (usually  signified 
by  broken   lines) 


X.         x^[ 


2.  Mark  m  points  z^Zp,...,z^  and  from  each  of  these 
draw  a  V  line  (usually  signified  by  a  broken  line 
to  a  small  filled  in  circle  marked  V) . 


V 


3.    Draw  directed  particle  lines,  one  entering  and  one 
leaving  each  of  the  points  x. ,  x-,  z..   These  lines 

X       -L       J 

may  run  between  points,  from  a  point  to  itself,  or 
from  a  point  to  the  edge  of  the  diagram  or  from  the 
edge  of  the  diagram  to  the  point. 

The  order  n  of  a  diagram  equals  the  sum  of  the 
interaction  lines  plus  the  V  lines,  i.e. 

n  =  p  +  m  . 

The  different  n-th  order  diagrams  are  obtained  by 
drawing  in  the  particle  lines  in  all  possible  ways. 

The  contribution  to  S  (t)  corresponding  to  a  given 
diagram  is! 
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1.  For  every  interaction  line  x.  xl    introduce  a  factor 
o-jj-  v(x.  -  X.  )  and  an  integration  with  respect  to 
dx.  dx!  over  the  region  of  space  time  behind  the 
surface  t. 

2.  For  every  V  line  whose  vertex  is  labelled  x. 

introduce  into  the  integrand  a  factor  -  — -  V(x.) 

and  an  integration  over  dx . . 

3 

3.  For  every  particle  line  from  a  point  y  (=  x.^x.  or  z.) 


to  the  edge  of  the  diagram  introduce  a  factor  ^"^(y) 


into  the  integrand.   For  every  particle  line  from  the 
edge  of  the  paper  to  a  point  y  introduce  f    (y)  into 
the  integrand.   If  f   (y)  corresponds  to  an  incoming 
line  and  f   (z)  to  the  corresponding  outgoing  line, 
then  Tp    (z)  is  adjacent  to  and  on  the  left  of  •^^(y). 

4.  For  every  particle  line  from  a  point  y  to  a  point  z 
(z  7^  y)  introduce  a  factor  S(z,y)  into  the  integrand. 

5.  For  every  particle  line  running  from  a  point  y  to 
itself  introduce  a  factor  p(y)  into  the  integrand. 

6.  Multiply  the  integral  by  — j-  (-1)^  where  n  is  the 
order  of  the  diagram  and  q  is  the  number  of  closed 
particle  loops. 

The  quantity  p(x)  is  the  charge  density  corresponding 
to  -pQ,    i.e. 

OCG 


1 
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The  u. (x)  were  defined  in  equation  (II. 6). 

u^[x)      =     u^(x)  e       .        -    ^ 

S(x',x)  is  a  propagator  which  is  a  solution  of*. 

ti^  5-|r  +k^'^  -  u(^')]  S(x',x)  =  5(5?-?')5(t-t-) 
and 


(2.9)     S(x',x)  =   9(t'-t)   ^_  u.(x')  a^(x) 

i 

occ 


-  9(t-t')  Yl'     u^(x')  a^(x) 


(t)  =  1   if  t  >  0 


0 

1 
2 


t  <  0 
t  -  0 


Note  that  whan  t  =  t',  S(x^x')  reduces  to  p-  -  p(x,x' ) 


where  p(x,x')  is  the  Dirac  density  matrix  corresponding 

to  the  state  f^,    i.e. 

occ 
p(x,x')   --=  Z^u^(x)  u^(x')  . 


To  clarify  the  above  prescriptions  consider  the 

o 


following  diagram; 


V 
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The  contribution  of  this  diagram-  to  S  (t)  is 


.  •^"''(x^)^"^(x^)  S  (x-]_,x^)  S(x2,x^)  S(x-^,X2)p(x^)  ]. 

4 .   The_  Linked  Cluster  Expansion . 

The  perturbation  series  for  the  energy  as  derived 
above  contains  terms  which  diverge  more  strongly  than  N^ 
the  number  of  particles  in  the  system^  as  N  — >  oo .  Such 
terms  can  have  no  physical  significance  and  it  should  be 
possible  to  eliminate  them  from  the  series.  This  is  the 
basis  for  the  linked-cluster  expansion. 

A  diagram  may  fall  into  two  or  more  unconnected  parts 
Such  a  diagram  is  described  as  unlinked.   Otherwise  it  is 
called  linked.   For  example! 


\ 


^- 


,->^^ 


unlinked  linked 
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It  is  the  unlinked  diagrams  which  have  no   physical 
significance  and  which  can  be  eliminated  from  consideration. 
It  can  be  shown  that 

S^(t)   ^   exp  {S^  (t)} 

a 
where 

linked 

a  i 

The  sum  is  over  all  linked  structures  F.   Note  that 

it  is  structures  which  are  summed  and  not  diagrams. 

Diagrams  which  have  the  same  arrangement  of  vertices, 

interaction,  V,  and  particle  lines  and  differ  only  in  the 

labeling  of  vertices  are  said  to  have  the  same  structure. 

S-p  (t)  can  be  resolved  into  two  parts! 
a 

(0) 
S^  (t)   =   S    (t)   +  S^  (t)  . 
a  a  a 

s|   (t)  contains  the  contributions  from  all  linked 
a 
vacuum  diagrams  i.e.  diagrams  without  particle  lines  going 

to  or  from  the  edge  of  the  diagram.   Hence  S-^   (t)  contains 

a 
no  operators  and  is  a  c  number. 

S^  (t)  contains  the  contributions  of  all  linked 
a 
diagrams  with  external  lines. 

It  follows  that*. 

1.  ip  =  Q.     lim  exp  (S^  (0))^q, 

c  is  a  normalization  factor. 

2.  whereas   AE  requires  only  a  vacuum  expectation; 
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Write  down  the  sum  of  all  integrals  contributing  to 

S^^  -^(00)   (one  from  each  vacuum  diagram).   Put  a  factor 

a 
iii5(t-j^)  into  each  of  these  integrands  and  put  a  =  0. 

The  resulting  sum  gives  AE. 


5  •   Polarization  Parts  _a_nd  _the_  Modified  Interaction . 
Consider  the  following  subdiagraml 


Diagram  A 

The  circle  represents  some  set  of  closed  particle 
loops  and  interaction  lines  which  are  connected  to  the 
rest  of  the  diagram  by  only  two  interaction  lines  as 
shown  in  the  diagram.   This  subdiagram  represents  the 
scattering  of  two  particles  through  some  complicated 
process  which,  however,  returns  the  particles  involved 
to  their  original  states.   Hubbard  shows  that  subdiagrams 
of  this  sort  represent  the  polarization  effect  and  can  be 
eliminated  from  the  perturbation  series  by  replacing  the 
ordinary  interaction  by  a  modified  interaction. 

The  term  polarization  part  is  used  to  describe  some 
connected  part  without  external  lines  which  is  connected 
to  the  rest  of  the  diagram  by  only  two  interaction  lines. 
For  example,  the  circle  in  the  above  diagram  is  a  polariza^ 
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tion  part. 

Suppose  r  is  some  structure.   If  F  contains  one  or 
more  polarization  parts  F  is  called  (polarization) 
reducible.   Otherwise,  F  is  irreducible.   All  reducible 
structures  can  be  obtained  by  inserting  polarization  parts 
in  place  of  interaction  lines  in  irreducible  structures. 
For  example,  the  above  diagram  A   can  be  obtained  from 
the  following  irreducible  one,  B,  by  replacing  the 
interaction  line  with  the  polarization  part  C. 


( 


V^' 


irreducible  polarization  part 

Diagram  B  Diagram  C 

Furthermore,  all  reducible  structures  are  uniquely 
reducible  i.e.,  any  reducible  structure  arises  from  a 
unique  irreducible  one. 

Consider  a  given  irreducible  structure  F'.  Then  it 
can  be  shown  that  the  total  contribution  to  S^i^t)  of  all 
diagrams  with  structures  which  reduce  to  F'  is  I 


S(t,F';v^) 
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where  this  expression  means  the  contribution  to  S(t) 

from  a  diagram  with  structure  T'  and  with  all  interactions 

V  replaced  everywhere  by  V,,  and 

(2.10)  V  (x',x)  ^      V::      W  (xSx.rJ) 

^         all  possible   ^       ^ 
polarization 
parts 

where   W, (x'^x^T  )  is  the  contribution  from  a  polarization 
part  Tg. 

Therefore  the  total  S-.(t)  is  a  sum  over  all  irreducible 
diagrams  with  v  replaced  by  V,  in  all  the  integrals.   Thus 
V,  playo  the  part  of  a  modified  interaction. 

For  calculating  the  wave  function  Vq  (i.e.  V  at  time 
equals  zero)  is  needed^  for  the  energy  shift  AE,  V    (time 
at  infinity)  is  needed. 

For  a  correct  counting  of  diagrams  one  gets  the 

following  expression  for  the  total  energy  shift  AE! 

irreducible  ^ 
va  cuum     „   ,  ■> 

(2.11)  AE  =       2_:        ^  i^(G)  AE  (G)  ; 

^      0 

n(G)  is  the  order  of  the  diagram  G. 

The  interaction  v  is  regarded  as  being  linearly 
proportional  to  some  coupling  constant  A. 

AE  (G  )  means  the  contribution  of  the  diagram  G  to  AE 
with  V  replaced  everywhere  by  V   .   (Henceforth  V   will 

r-  '^  <^    CO      ^  OD 

be  designated  as  V. ) 

The  sum  is  over  all  irreducible  vacuum  diagrams. 
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6  .   AiT^_Int_e_gra_l  Equation  for  _th  e  Modified  Interaction . 

The  modified  interaction  V  has  been  defined  as  an 
infinite  sum  over  all  polarization  parts.   This  is  not 
a  rapidly  convergent  series .   It  will  now  be  shown  that 
there  exists  an  integral  equation  for  V  and  that  this  allows 
us  to  calculate  V  from  a  series  which  is  rapidly  convergent. 

A  polarization  structure  which  consists  of  two  or 
more  parts  which  are  connected  only  by  single  interaction 
lines  is  called  an  improper  polarization  part.   If  the 
structure  cannot  be  resolved  into  at  least  two  parts 
which  are  connected  by  a  single  interaction  line,  it  is 
called  a  proper  polarization  part.   For  example,  diagram  D 
is  a  proper  polarization  part  and  diagram  E  is  an  improper 
one. 


X 


o 


D  S 

Suppose  r  is  some  improper  polarization  part.   Then 
it  can  be  uniquely  resolved  into  a  proper  part  and  some 
other  polarization  part,  proper  or  improper. 

Let  W(x',x,r)  stand  for  the  contribution  to  S  (t) 
from  the  diagram  F.   Define  ¥(x',x  T)  by! 
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2.12)   W(x',x,r)  =   v(x'-x")  W(x",x"Sr)  v(x"'-x)  dx"  dx'" 


The  two  v's  correspond  to  the  two  outgoing  interaction 
lines  of  the  polarization  structure  and  W  corresponds  to 
the  rest  of  the  structure.   See  diagram  F  as  an  example. 


V 


X 


--0-- 

X  "  ^'^  X  "  ' 


V 


X 


Diagram  F 

The  whole  diagram  contributes  W  and  the  part  within 
the  circle  contributes  W.   For  the  example  shown  in 
diagram  D,  W  is 


(2.15) 


W  =  -  ^-  S(x',x)  S(x,x') 


ih 


If  the  improper  structure  F  is  resolved  into  a  proper 
part  F'  and  another  part  F ", 


(2.14)   W(xSx,F)  .  J  v(x'-x^)  W(Xj,X2.r')  ^(x^-x^) 


•  W(x^,X|^,r")  v(x^-x)  dx^  dx2  dx^  dX|^, 


Define  w'^(x',x,r')  as 


(2.15)   ¥'''(x',x,r')  =  f  v(x'-x")  W(x",x,F')  dx" 


and  (2.l4)  becomes 


r 


(2.16)   ¥(xSx,r)   =   !  W*(xSx",r')  W(x%x,  F")  dx" 
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The  total  contribution  to  V  from  all  improper  structures 
which  resolve  into  the  proper  part  F'   is  obtained  by 
summing  over  all  structures  T ".   From  the  right  side  of 
(2.16)  this  sum  is 


(2.17)   fw^  (x-.x-'^rt)  —   (^-'^^^  r„)  ^^n 

F" 

s=    r  W*(x',x",F')  [v(x"-x')  +  X~  W(x",x,  F")]  dx" 

J  riTTn 


F'VF 

where  v(x"-x')  is  equivalent  to  F    But  the  sum  over  all 
polarization  structures  F"  is  just  the  definition  of  V. 
Thus  (2.17)  can  be  written  as 

(2.18)        j  ^*{x^,x\r^)   V(x",x)  dx"  . 

This  is  the  contribution  to  V  of  the  proper  part  F' 
and  all  the  improper  parts  which  resolve  into  F '.   To  obtain 
the  total  V  ^  (2.l8)  must  be  summed  over  all  proper  polariza- 
tion parts  F'. 

Let 


proper 

(2.19)  v*(xsx)  -     y       \^*{x^,x,  r-) 

T' /Fq 


so 


(2.20)   V(x',x)   -  v(x'-x)  +  f  V*(xSx")  V(x%x)  dx"  . 
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7 .   Tti_e  Self-Gonoistent  Field . 

The  remarks  in  this  section  apply  only  to  vacuum 
diagrams  (those  without  external  lines).   As  has  been 
discussed  previously^  these  are  the  only  diagrams  which 
contribute  to  the  energy  shift. 

Any  part  of  a  diagram  which  is  connected  to  the  rest 
of  the  diagram  by  only  a  single  interaction  line  is  called 
an  H-part.   For  example,  diagrams  G  and  H  are  both  H  parts. 

k — xO         '"'O^ "h3 

G  H 

If  there  is  no  external  potential  then  it  can  be  shown 
that  diagrams  containing  H  parts  do  not  contribute  to  the 
energy.   The  proof  depends  upon  the  uniformity  of  the  gas 
and  so  does  not  apply  when  there  is  an  external  potential. 
However^  the  potential  V  can  be  chosen  in  such  a  way  so 
that  the  H  parts  can  be  eliminated. 

If  the  H  part  can  be  resolved  into  a  polarization  part 
and,  some  other  H  part,  it  is  called  improper.   Otherwise 
it  is  proper.   For  example  diagram  G  is  proper  and 
diagram  H  is  improper. 

An  H  part  F  gives  a  factor  which  we  will  denote  by 
H(x.,r).   There  is  no  preferred  origin  of  time  in  the 
theory  so  that  H(x.,r)  is  independent  of  time,  and  can  be 
written  as  H(x.,r). 
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We  can  choose  V(x)  so  that  the  subdiagram  I  cancels 
the  sum  of  all  subdiagrams  J.   (in  diagram  J,  F  is  a 
proper  H  part . ) 


)v-© 


If  a  correct  counting  of  diagrams  is  performed,  it 

can  he  shown  that  the  choice  for  V(x}  is 

on(r) 
(2.21)         V(5?)  =  ^::  2__H'(x,r)  . 

The  sum  is  over  all  proper  H  parts*  n(r)  is  the  order 
of  the  H  pari:  r,  g(r)  is  the  number  of  permutations  of  the 
interaction  lines  which  leave  F  unchanged. 

With  this  choice  for  V  one  can  eliminate  all  diagrams 
containing  H  parts  with  the  exception  of  symmetrical 
diagrams  of  the  type  shown  in  diagram  K. 


0-  •-0 


K 

The  total  contribution  to  the  energy  from  these  types 
of  diagrams  can  be  shown  to  be 
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_  p2n(r)-l 
(2.22)   Ey  =  -  X_  ^ ^J  H(x,r)  v(x-x')  H(x',r)  dx  dx  •  . 

r  [g(^)] 

The  simplest  H  part  is  that  shown  in  diagram  G.   If  all 
other  H  parts  are  neglected,  we  have 

(2.25)         V(x)  =  Jv(x-x')  p(x''')  dx' 

which  is  the  same  as  the  symmetrized  Hartree  field. 
Therefore,  the  self-consistent  field  introduced  in  this 
section  is  a  generalization  of  the  ordinary  Hartree  field. 
Keeping  only  the  simplest  H  part, 

(2.24)    Ej^  =  _|.Jp(^)  v(x-x')  p(x')  dx  dx'  , 

which  is  the  term  which  must  be  subtracted  when  calculating 
the  total  energy  in  an  ordinary  Hartree  calculation. 

The  calculation  of  V(x)  from  (2.21) is  a  self-consistent 
calculation  since  H'(x,r)  depends  upon  the  density  p(x) 
and  the  propagator  S(x',x).   These  depend  upon  the  u. (x) 
and  E.  which  in  turn  depend  upon  V(x).   (See  equation  (2.6).) 
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8 .   The  Energy  Shift  and  Collective  Modes . 

Making  use  of  equatioa  (2.11)  and  of  correct  counting 
procedures,  it  can  be  shown  that  the  energy  shift  AE  is 

(2.25)  AE  =  i^  J  ^  J  dx^  J  dx{  6(t^)  V(x^,x|)V(x;,x^) 

0 

—         proper^ 

(2.26)  V(x',x)  =   ^__  W(xSx  r) 

where  W  has  been  defined  in  equation  (2.12)  and  the 
definition  of  proper  is  discussed  just  prior  to  equation 
(2.12). 

From  equations  (2.I5)  and  (2.19)  we  have 

(2.27)  V*(x',x)  =  Jv(x',x")  v(x"-x)  dx 
and  (2.20)  can  be  used  to  calculate  V(x«,x-,) 

(2.20)   V(x'  x-^)  =  v(x|-x^)  +  J  V*(x|,x")  V(x%x^)  dx"  . 

The  total  energy  of  the  system  is 

A 

OGG 

(2.28)  E  =  y;^  E^  -  i-  Nv(0)  +Ej^+^J^Jdx  dx' 

^  0 

•  5(t)  V(x,x')  V(x',x)  . 

In  order  to  study  AE  we  must  solve  the  integral 
equation  (2.20).   For  a  uniform  gas,  this  solution  is 
easily  done  by  Fourier  transformation.   For  the  non- 
uniform case  we  need  a  different  method  which  will  now 
be  described. 
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Multii-  ^   ,       ,  ,  ,      ^     .,x^x")  and  integrate 
over  x".   Than 

{2.2Q)      F(x',x)  =  V*(x',x)  +  rv*(x'x,")  F(x".^  '  •: " 


(2.?0)     F(x',x)   =   J  V(x',x")  V(x".  ^  -." 

I,  .  ,)    V*(x],x)  =  Jv(x',x")  v(x"-x)  dx" 

.  .   )     AE  =  ^  J"  ^Jdx  6(t)  F(x,x)  . 

0 

It  is  still  possible  to  Fourier  transform  with 

respect  to  the  time  difference. 

(2.51)  F(x',x^cD')  r--  V*(x«,x  00')  +  J  V*(x',x",co')F(x",x  a3')dx" 

(2.52)  AE  =  j^  J^  /do.-  dx"F(3^,x^co.) 

0 

To  solve  equation  (2.29),  introduce  the  eigenvalues 
and  eigenfxinctions  of  the  homogeneous  equations 

Hj_(cu')4>^(x,aj' )  =  J  V*(x,x',a)')l>^(x',OD' )  dx' 

(2.33) 

^^i(<J^')X(x,co' )   =  J  V*(x',x,cu' )x^(xSco')  dx'  . 
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If  F(x',Xja)')  is  expanded  in  terms  of  the  4>.(x')  and 
6  "■ 

Tp..{x)  and  use  is  made  of  the  orthonormality  relations 
J 

between  the  ^.(x)  and  <i> .  (x)  ,      we  can  solve  equation  (2.51) 
and  equation  (2,32)  reduces  to 


A  /   h.  (oj'  )  \i.  (go'  ) 

n  .3  ^  a  J      -L  J- 


^^      ^   do)'   i 


^0  1  -  '^("^') 

where  h.  (00' )  is  the  degeneracy  of  the  eigenvalue  ix.  (cd' ) . 
Writing  1-l.  (cd' )  =  A.(cd')  +  i  2.  (cd' )    and  remembering 
that  AE  must  be  real 

t.  ?dA  r       ^i  ^(^')  ^i^"^') 


I     ""      ^    [1-A.(CD.)]-  -f  2^(0)') 


5""": 


Since  V*  (x»  ,x",cd' )  is  not  necessarily  a  self -adjoint 
operator^  we  realize  that  the  solutions  <t>.  and  X^  of  tne 
homogeneous  equations  may  not  form  a  complete  set. 
However,  we  note  that  any  operator  can  be  put  into  Jordan 
canonical  form  and  also  that  (2.31)  can  be  written  as 
F  =  (1-V*)"-^V^.   In  (2.32)  we  are  interested  in  Tr  F. 
Using  these  facts  we  can  show  that  the  lack  of  a  complete 
set  simply  causes  a  modification  of  the  weighting  factor 
h(oo')  in  (2.3^). 
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There  will  iDe  a  large  contrilDution  to  the  energy  at  those 
irequencies  co'  for  which  any  A.(a)')  =  l  provided   ^{^' ) 
is  small.   These  large  contributions  to  the  energy  can  be 
associated  with  collective  modes.   The  frequency  of  these 
modes  is  given  by  A.(co  )  =  1  and   ^{^^)    describes  the 
damping.   If   2(co  )  is  large  it  will  not  be  possible  to 
recognize  the  collective  modes.   Thus  the  conditions  for 
a  collective  mode  are 

A^(a3')   .   1 

(2.56) 

2.  (CO-)  .   0  . 

We  will  show  more  explicitly  in  the  next  chapter  that 
equations  (2.56)  really  do  correspond  to  a  collective  state 
Also,  in  Chapter  IV  we  will  see  that  these  equations  lead 
to  an  accurate  prediction  of  the  collective  modes  for  a 
free  gas. 


-27- 


CHAPTER  III 

7 

Physical  Interpretation  of  the  Collective  Mode ' 

1 .      General   Rema rks  . 

Experimentally,  one  measures  the  excitations  of  a 
many  body  system  by  means  of  an  external  probe  such  as 
inelastically  scattered  particles  or  X-rays.   For  example, 
the  experiment  might  measure  the  energy  distribution  or  the 
angular  distribution  of  the  scattering  cross-section. 

If  the  density  of  the  many  body  system  is  p(r)  we 
can  write 

/  \   V A/     \    '^: ik-r    -r ^i  k  •  (r-r.  ) 

p(r)  =  > 6(r-r.  )  =  L_  Pv  ^      =  ^__  ^    ^   i"' 

i       ^     k   ^  i  k 

where  p,  is  called  the  density  fluctuation  and  is  defined  as 

n:^ —  ^i  k'  ri 

Pk  =  2__  e 

This  p.  describes  the  fluctuations  in  particle  density 
about  the  average  density.   The  probability  that  the 
external  probe  transfer  momentum  k  and  energy  o)  to  the 
many  body  system  is  proportional  to 


>    (pM^n  <5(ai-a)  ^  ) 
~ ^^  k  'nO   ^   nO  ' 


where    the   co  ^  are   the   exact   excitation  frequencies   of   the 


'    General   references    for   this   chapter  are  I 

D.    Pines,    The   Many   Body   Problem,    pp.    I-69,    New  York    (1962) 
D.    Pines,    Elementary   Kxcitations    in  Solids,    pp.    55-155 
New   York,  (196^}". 
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many  body  system  induced  by  the  density  fluctuations  pj^. 
The  proportionality  factor  depends  only  on  the  interaction 
v(r)  between  the  incoming  probe  and  the  target  particle 
and  not  on  any  dynamical  properties  of  the  many  body  system, 

For  example,  if  the  Born  approximation  applies,  this  factor 

2 
would   be    27rv,    where 

r      ,    N      ik-r    ,3 
v^     =     J    v(r)    e  d^r    . 

The  maximum  experimental  information  about  the  many 

body  system  is  therefore  contained  in  the  quantity 

2 
S(ka3)  =  5    |p+|  _  6(00-0)  \    . 

S(k(jL))  is  the  Fourier  transform  in  space  and  time  of 
the  time  dependent  pair-distribution  function  for  the  many 
body  system.   It  is  the  spectrum  of  excitations  available 
to  the  density  fluctuations  p   that  is  actually  measured, 
i.e.  the  excited  states  are  coupled  to  the  ground  state 
by  the  density  fluctuation  p,  . 

For  a  non-interacting  fermion  system,  the  density 
fluctuation  spectrum  coincides  with  simple  single  particle- 
hole  pairs. 

Consider  now  a  fermion  system  with  an  interacting 
potential  between  all  the  particles.   This  interaction 
causes  correlations  in  particle  positions  and  these 
correlations  cause  coherence  effects.   In  a  dense  electron 
gas,  the  correlations  give  rise  to  screening  i.e.  each 
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electron  polarizes  its  immediate  surroundings.   The  density 
fluctuation  excitation  spectrum  tajIII  contain  instead  of 
simple  particle  hole  pairs,  quasi-particle  hole  pairs. 
The  quasi-particle  is  the  particle  plus  a  screening  cloud. 
In  addition,  sometimes  it  is  possible  that  the  correlations 
give  rise  to  coherence  factors  which  greatly  reduce  the 
probability  for  exciting  single  quasi-particle  pairs. 
Instead,  there  may  be  a  coherent  superposition  of  particle 
hole  pairs  which  is  referred  to  as  a  collective  mode. 
The  collective  state  differs  from  the  quasi-particle  states 
both  in  amplitude  and  spacing.   The  amplitude  of  a  collective 
state  is  proportional  to  N,  the  number  of  particles  in  the 
system.   Single  particle  amplitudes  are  not  proportional 
to  N.   Also,  for  a  given  momentum  transfer  we  would  expect 
collective  states  to  be  widely  spaced  whereas  single  particle 
states  are  close  together.   Also,  since  it  is  a  coherence 
of  density  fluctuation  excitations,  the  collective  mode 
corresponds  to  oscillations  in  the  particle  density.   For 
this  reason,  when  speaking  of  free  gases,  collective  modes 
are  sometimes  referred  to  in  the  literature  as  zero  sound. 
However,  unlike  sound  waves  it  is  not  the  frequent  collision 
of  particles  which  acts  as  the  restoring  force.   The  restoring 
force  is  the  averaged  force  field  of  a  large  niunber  of 
interacting  particles.   Collisions  between  the  particles 
in  a  gas  would  disrupt  this  averaged  field  and  damp  the 
collective  mode.   Therefore,  we  would  expect  collective  modes 
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to  be  observed  only  if  cut  >>  i,    where-  T  is  the  average  time 
between  aollisions.   This  condition  is  satisfied  at  low 
temperatures.   In  this  paper,  we  assume  the  temperature  to 
be  zero,  so  the  above  condition  is  always  satisfied.   This 
assumption  is  valid  for  a  system  such  as  an  atom  at  room 
temperature.   However,  in  a  stellar  system  where  the  tempera- 
ture is  very  high  the  results  of  this  paper  would  not  be 
valid.   Furthermore,  we  might  expect  that  in  a  star  we 
would  not  observe  any  collective  modes  as  they  might  be 
heavily  damped  due  to  particle  collisions. 

For  collective  modes  to  exist  and  to  be  observable 
there  is  another  requirement.   They  must  not  decay  quickly 
into  single  quasi-particle  pairs.   Either  this  decay  must 
not  be  energetically  possible  or   else,  because  of  coherence 
effects  between  the  relevant  matrix  elements,  the  probability 
for  exciting  these  single  particle-hole  pairs  must  be  greatly 
reduced.   If  this  probability  is  high,  the  collective  mode 
would  be  very  quickly  damped  and  so  'on observable.   If  this 
probability  is  zero,  then  the  collective  mode  can  be  a  long 
lived  well  defined  excitation  of  the  system.   The  damping 
term  referred  to  in  this  paper  is  a  measure  of  this  proba- 
bility for  exciting  single  quasi-particle  hole  pairs. 


-31- 


2 .   Relation  to  the  Dielectric  Constant. 

The  Fourier  transforms  In  time  of  the  macroscopic 
Polsson  equations  are: 

div  D(a))   =   £(co)  div  E(a))  =  47rq(a)) 
div  E(tD)   =   4(r[q(cD)  -  e<p(a))>] 

where  q(a))  is  the  true  free  charge  and   e<^p(co)">  represents 
th:  polarization  charge.   If  q(<^)  represents  an  external 
prohe,  then  p(<Jo)  represents  the  response  of  the  system. 
A  collective  mode  corresponds  to  a  free  resonance  in  the 
system,  i.e.,  in  the  absence  of  any  external  q  there  must 
he  a  fluctuating  non-zero  electric  field  E(co)  and  also 
a  non-zero  fluctuation  in  \p(<^),'>.   The  only  way  for  E  and  p 
to  be  nonzero  if  q  is  zero  is  for  some  eigenvalue  of  the 
dielectric  constant  operator  to  be  zero.   So  a  collective 
mode  corresponds  to  the  condition  that  e(<^)  is  singular. 

o 

It  has  been  shown  by  Nozieres  and  Pines   that  the 
shift  in  the  total  energy  of  the  system  due  to  the  inter- 
action between  the  electrons  can  be  written  in  the  form'. 

oo  2 

0 

where  by  Im  we  mean  the  skew  hermitian  part. 


P.  Nozieres  and  D.  Pines,  II  Nuovo  Cimento  [x],  %    (1958), 
pp.  ^70-^B9. 


-32- 


Actually  Nozieres  and  Pines  derived  this  relationship 
for  a  uniform  system  where  e(a))  is  diagonal  in  a  plane  wave 
representation.   We  believe  their  proof  can  be  generalized 
so  that  a  relation  very  similar  to  (j).l)  is  true  for  a 
non-uniform  system. 

We  must  subtract  the  term  Ne  /r  so  that  this  definition 
of  the  energy  shift  coincides  with  the  energy  shift  defined 
in  equation  (2.25).   From  equation  (2.34)  we  have 

(2.54)  AE  =  Re  [^  ^  ,  f  J  d..  _%i_2^L^ 

^^;here  h  (  oj  )  is  the  degeneracy  of  each  M- ( oj  '  ) . 
It  is  not  hard  to  show  from  (3.1)  that 

0 

By  comparing  these  two  equations ^  we  see  that  the 
condition  for  a  collective  mode  that  was  previously  stated, 
i.e.  that  M-(a3' )  =  1,  is  completely  equivalent  to  the  condi- 
tion that  £(a))  be  singular. 

We  think  that  a  collective  state  in  a  non-uniform  system  is 
a  new   type  of  excitation  and  as  such  would  not  show  up  by 
perturbational  treatment.   This  means  on  the  one  hand  that 
if  the  non-uniformity  is  treated  perturbationally  with 
respect  to  the  free  gas  the  result  will  be  the  same  as  in 
the  free  gas.   We  show  this  in  Appendix  C.   On  the  other  hand. 
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Atlick  and  Glassgold   have  shown  that  the  correlation 
energies  for  single  particle  states  when  calculated  in. the 
random  phase  approximation  are  quite  small  corrections  to 
the  Hartree-Fock  energies.   Therefore,  we  think  that  a  very- 
large  contribution  to  the  correlation  energy  must  correspond 
to  a  new  type  of  excitation,  namely  a  collective  state. 

5 •   Structure  Factor  _S  ( k )  and  the  Spectral  Function  S  ( k , oo ) 

Pines  has  derived  a  useful  relationship 

P  p 

^^  TTklFT  =  -  —J?-  s(i^^)  ^    0.  >  0  . 

Thus  a  zero  in  e(kcu)  should  correspond  to  a  peak  in 
S(kco).   (S(kcD)  has  already  been  defined.) 

The  structure  factor,  S(k),  =  ^^^o'^k  ^kl  ^(/^  /  ^ 
and  is  the  Fourier  transform  of  the  time-independent  pair 
distribution  function. 

S(ic)  is  related  to  S(lca))  by 

00 

^  dco  S(foi))   =  S(k)N  . 

V 

-00 

If  we  refer  back  to  the  introduction  we  see  that  there 
may  be  more  than  one  discrete  eigenvalue  |x.  (cd).   Each  of 
these  eigenvalues  may  yield  one  or  more  (or  perhaps  none) 


^Atlick  and  Glassgold,  Phys .  Rev.  l^l  (1964)  A652-A646. 
■'"'"'d.  Pines,  The  Many  Body  Problem,  p.  ^Q,    N.  Y.  (I962). 
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collective  modes.   If  we  knew  the  complete  solution^  i.e. 

all  the  collective  modes  and  all  the  eigenvalues  we  would 

have  the  equivalent  of  a  discrete  dispersion  relation. 

In  the  free  gas  case  the  index  1  becomes  equivalent  to  the 

continuous  wave  vector  k  ,  and  M-.  (cu)  is  replaced  by  H.(kcL)). 

The  collective  mode  condition  1-  (kco)  =  0  provides  a 

dispersion  relation  <^(k).   Wa  do  not  have  such  complete 

information  here  as  we  have  evaluated  only  one  eigenvalue. 

We  know  at  what  frequency  co  to  expect  a  collective  resonance 

but  we  do  not  have  a  dispersion  relation.   Therefore,  we 

know  that  S(kci^)  has  a  peak  at  Oi     but  we  do  not  know  the 

complete  form  for  S(ka))  and  so  we  cannot  evaluate  S(k). 

These  remarks  can  be  illustrated  more  clearly.   R.  D.  Puff  [15] 

has  shown  that  for  a  single  resonance  <^Q(k)  with  zero  damping 

There  is  a  sum  rule 

OD 

I  J  g-  coth  I  pa3(2|)  s(ka3)   =  ^^  ,   where  ^   =  ^^  , 

-00  ' 

With  the  above  S(ka3)  this  sum  rule  becomes! 

2 

If  p-  po)  (k)  >>  1,    which  is  certainly  true  for  the 
temperature  and  energy  ranges  we  are  working  with,  this 
S'om  rule  becomes  „ 


S(k^)  =  ^fk-f"  ^^(^-^0^^))  +  ^("^  +  ^0^^))  • 
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Thus  tha  structure  factor  is  determined  by  the 
dispersion  relation  which  wa  do  not  taow.   Actually,  we 
could  calculate  S(k).   In  his  paper.  Puff  derives  a  sum 
rule  which  provides  an  integral  equation  for  S(k). 
However,  as  will  be  shown  in  the  next  section,  S(k)  is 
not  needad  for  comparison  with  experiment  as  S(ka)) 
suffices.   Therefore,  we  will  not  try  to  calculate  S(k). 

4 .   Relation  to  Cross  Sections. 

If  particles  are  scattered  from  the  many  body  system 
then  the  scattering  cross  section  is  given  by 

d^'-?c5  -     ^     ^-kS(^-)      (van  Hove)  , 

where  it  is  assumed  that  the  Born  approximation  applies. 

M  is  tha  mass  of  tha  particle  being  scattered^  Fj   is  tha 

initial  momentum,  P^  the  final  momentum,  with  energy 

transfer  between  cd  and  cd  +  dcu  and  sca"t"'"ering  angle  dO. 

ik*  r 
V^  =  J  e     V(r)  d^r  and  ii  =  1. 

If  we   integrate    over  all   energy   transfers 

00  2      p 

dQ      -      J    "^"^  dco  dO      -      o    5      P^    \  ^^^^^    • 

Thus  we  see  that  if  the  scattering  cross  section  is 
measured  as  a  function  of  the  energy  transfer  co  for  a 
fixed  scattering  angle,  the  cross  section  should  have  a 
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peak  at  those  co  corresponding  to  the  qollective  mode. 
If  the  angular  distribution  of  the  scattered  particles 
is  measured  without  specifying  energy  transfer,  then  S(k) 
is  determined. 

We  should  mention  that  since  in  this  work  we  required 
zero  damping  the  peak  in  S(kco)  is  infinitely  high.   If  the 
problem  were  done  more  carefully,  for  example,  if  more 
diagrams  were  included  in  the  evaluation  of  V,  the  damping 
would  not  be  zero  and  the  peak  would  have  finite  height. 

Normally  if  X-rays  are  scattered  from  the  system,  then 
the  observed  cross  section  is 

(5.3)  I?  =   ao-s(K)  , 

where  d  a^  /  dO   is  the  differential  scattering  cross  section 

2 
for  X-rays  and  a  single  electron.   At  energies  <<  mc  , 

where  m  is  the  mass  of  an  electron,  da^  /dO  reduces  to  the 

Thomson  cross  section  d  a„  /dO  . 

In  usual  X-ray  scattering  experiments,  the  energy 

transfer  is  small  compared  with  the  experimental  resolution, 

and  the  observed  cross  section  does  not  depend  on  the  energy 

transfer.   However,  since  we  know  more  about  S(kaj)  than 

about  S(k),  it  is  more  convenient  to  look  at  the  differential 

cross  section  with  respect  to  the  energy.   Since  the 

derivation  of  the  Thomson  scattering  cross  section  is 

independent  of  the  frequency  we  have 


-37- 


dO  d'X)     dO    ^  ' 


^      i  (-%)  (2  -  sin2  e)  . 


The  cross  section  should  show  a  peak  at  the  frequency 
of  the  scattering  X-ray  which  corresponds  to  o^^,  the 
predicted  collective  frequency.   Experimentally,  this 
might  be  determined  from  the  Mossbauer  effect. 
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CHAPTER  IV 
Hubbard's  Solution  for  a  Uniform  Gas 


There  are  some  important  reasons  for  including  the 
solution  for  a  uniform  gas  in  our  discussion.   The  solution 
for  the  uniform  gas  guided  the  approximations  used  in  the 
solution  of  the  more  general  case.   Furthermore,  the  general 
solution  must  reduce  to  the  free  gas  solution  when  the 
external  field  goes  to  zero.   Also,  it  is  interesting  to 
see  how  much  simpler  the  problem  is  when  there  is  no 
external  field. 

Consider  now  a  many  body  system  of  electrons  with  no 
external  forces  present.   The  interaction  v(x-x' )  is  the 
Coulomb  interaction.   In  Chapter  II  it  was  found  that 

A 


(2.25)  '^E  =  ^  /  f-  /  dx^  J  i^l    6(t^)  V{.^,^[)    V(x;,xj 


0 
where 


(2.20)   V(x^,x^)  =  v{x[-x^)    +   J  V*(x^,x")  Y{x'\x-^)    dx" 
v(x^-x^)  =  v(x^-x^)  6(t-J^-t-^) 


and 


^    rr/ 


(2.27)    V  (x',x)=  j  V(x',x")  v(x"-x)  dx"  , 
and 


11 


J.  Hubbard,  Proc.  of  the  Royal  Soc.  (London)  A243 
(1958)  336-351.  — 
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_  prop9r_ 

(2.26)         V(xSx)   -   '^~  W(x',xr)  . 

In  the  absence  of  any  external  fields  V_,  V*  and  V 
are  functions  of  only  the  differences  in  their  arguments. 
Therefore  equation  (2.20)  can  be  solved  immediately  by 
introducing  the  Fourier  transforms  of  V^  V   and  Y', 

-> 

(4.1)  v(xSx)=i:-  J  dk  J  dc  V(k,co)  ei^-(^'-'^)+i^(^'-*) 

The  definitions  for  V  (k^co)  and  V(k,co)  are  similar. 
Equations  (2.20)  and  (2.27)  become 

— > 

(4.2)  V(kjco)  =  v(l?)  +  V*{^jCo)    7(1^^03) 

(4.5)  V*(k/D)  =  v(k)  V(k/o) 


i^A)  v(k)   =  ^f-. 


2 
k^ 

which  is  the  Fourier  transform  of  the  Coulomb  potential. 
Equation  (4.2)  can  be  solved  immediately. 


(4.5)  v(k,6o)  = 


vM 


1  -  V*(k,a)) 
The  expression  for  AE  becomes 


(4.6)         AE-=i^-   f^  r^k     P 


47rn  J  A  J  ^27r)^  ^    1  -  V*(k,co) 


where  we  have  divided  by  the  number  of  particles  in  the 
system  and  n  is  the  particle  density. 
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^      — >  — >  — f 

Writing  V  (k,oo)  :=  A(kja))  +  i  51  i^}^)    > 


^  ■  ^   (k,^) 


0 


^Fn  J  ^  J  (2Tr)^  ^     [1-A(l?,a3)]^+  ^^(i?,co) 


From  (2.26)  we  have 

_  -»      proper  _  -» 
(4.8)  V(k,a))  =  ^   W(k,a)) 

where  the  sum  is  over  all  proper  polarization  parts. 
This  sum  is  approximated  by  its  leading  term,  that  corres- 
ponding to  the  simplest  of  all  polarization  parts  shown  In 
the  following  diagram. 


The  expression  for  W  corresponding  to  this  diagram  is 
(4.9)  w"(xSx)   =.   -  ^  S(x',x)  S(x,x')  . 

For  a  free  gas 


unocc  (ik*  (x'-x)-  J^  E,  (t'-t)) 
(4.10)   S(x',x)   =  9(t'-t)  ;^  e  ^     ^ 

k 


occ   (ik- (x'-x)-  i  E.  (t'-t)) 
-9(t-t')   ^e  ^  ^ 


where  E^  =   ii^k^/2m. 
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SulDstitating    (4.10)    into    (4.9)   we   get 


mocc  ocG 
(4.11)    W(x'-x)    =  ^   [9(t'-t)   TZlZ 


i(k'-k")-(x'-x)-   i(E^,-E^„)(f-t; 


oc_c   unqcc      [i  (k' -k") .  (x« -x)  -  ^(E,  , -E,  „  )  (t ' -t)  ] 

+   9(t-t')    N 1_        e  n      K        IV  ^ 

k'         k'^ 


¥(x'-x)  is  Fourier  transformed  to  give  W(k^cD). 

■X-  —  "^         "^   47re^ 

Then  V  (k,(X3)  is  derived  by  multiplying  W(kja3)  by  v(k)  = 

k' 

The  result  is 


_e_ 
2 


r^/T^  ..\   4Tr  e^ 


(4.12)  V  (k.cu)  =  ^"-"-  >__  N -2 r2Pr- 

^    m      '^      4m^ 


^'•-'^-¥  >=  \.(1  -  ''...K)tM'^'^K---  -'( 


ii  k    k' 


/ 


/'l  if  state  k  is  occupied 
where  P  means  principal  part  and  N,  =  ( 

^   lo  "   "     "  not  " 

The  real  part  of  this  expression  is  A  and  the  imaginary 

part  is  2  . 
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Letting  x  •=  ;     .nd  y  -  h  cd /cl    it  is  found  that 
r  r 

for  small  x  and  for  y  >>  x(x+2) 

2(x,y)   =   0 

2 
A(x,y)   =  -^ 

y 


V 


Thus  A(x,y)  =  1  irr;  '  "  -.  y  =  y  and  the  frequency  equals  the 

It 

plasma  fr^     :y. 

Be      leaving  this  .;     a,    it  is  interestii._ 
compare  equation  (^.6)  for  t!     'rgy  shift  for  a  uniform 

■3tion  (2.34)  which  is  the  energy  shift  for 
a  general  system. 

{2.:  \  r  ¥  f  d^'  — — 

^'   i^-^  1  -  ^x.(a).) 

where  the  M-.  are  the  eigenvalues  of  the  homogeneous  integral 

;  .-.ition 

~*        p  *  ~*  -^ 

M.^(oD')  <l>^(x,oo')  =  J  T(x,x',a3')  «J)^(x',u)')  dx'  . 

»Ve  see  that  for  a  uniform  gas,  the  eigenvalue  |i.=  M-Ck) 
is  just  the  Fourier  transform  of  V  (x^x'^od'),  and  the 
eigenfunctions  <l>.(Xja)')  =  <{>(ic,x,tt)» )  are  just  plane  waves 
or  weighted  spherical  averages,  i.e.  free  spherical  waves. 
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CHAPTER  V 
^^?-  Solution  for  a  Dense  Electron  Gas  in  an 
External  Potential  -  General  Remarks 

It  was  shown  in  Chapter  II  that  the  information  about 
collective  modes  is  contained  in  the  eigenvalues  M-  (co)  of 
the  integral  equations 

(2.33)    ^i^(a3)*^(x,co)    =   J  V*(x,x',a3)(D^(x',a3)    dx ' 

h  ("^)X,  (x,co)    =    f  V*(xSx,co)X.  (x',cd)    dx' 
where 

(2.27)    V*(x',x)      -     Jv(x',x")    v(x"-x)    dx" 
/rr   1  \  ^    II      \        5(t"-t)e^ 

(5.1)  v(x"-x)  =  -y—f, — 


TT — n 
X    -x 


and  V  was  defined  as  the  sum  of  proper  polarization  parts, 
(See  Chapter  II,  Sections  5  and  6).   For  the  calculation 
of  V,  we  shall  keep  only  the  leading  term  -  the  term  that 
corresponds  to  the  simplest  of  polarization  parts! 


It  can  be  shown  that  V  is  a  rapidly  converging  series 
for  a  weak  interaction  v  or  a  very  dense  gas.   Therefore _, 
the  keeping  of  only  the  leading  term  should  be  a  good 
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approximation  for  a  high  density  ele.ctron  gas.   This 
approximation  is  often  referred  to  in  the  literature  as 
the  random  phase  approximation;  sometimes  it  is  called 
the  pair  approximation.   We  saw  that  it  was  quite  successful 
in  predicting  the  collective  motion  for  a  free  gas. 
Therefore 

(5.2)  V(xSx)  =  -  3^^  S(x',x)  S(x,x') 

where 

unocc 

(5.3)  S(xSx)  =  9(t'-t)   ZZ  ^k(^')  u  (x) 

-  0(t-t')  Yr_  u^(x')  ^^^(x) 

k 
where  9   is  the  step  function  and  the  u,  (x)  were  defined 
in  equation  (2.6). 

(2.6)   [£^  +  U(^)  +  V(^)]u^(3?)   =   E.u^^(^)  . 

U(x)  is  the  external  potential  and  V(x)  is  the  self- 
consistent  field.   In  the  random  phase  approximation  V(x) 
is  the  Hartree  field.   Henceforth,  we  will  assume  V(x) 
known  and  refer  to  the  sum  of  U(x)  and  V(x)  as  U(x). 
From  (5.2)  and  (5.3)  we  have 

unocG  occ 

k   '^      '^     i 

•0(t'-t) 


—         -,    unocG  occ 

(5.4)   V(x',x)  =  ^^-  [(  ^  Uj^(x')  n^ix)   JZ   ^i(x)  u.(x-)) 


unocc  occ 

+  e(t-t-)  XI  u.  (x){I.  (x')^  u^^(x')  u^{x))] 


k 
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Now  J 

(5.5)   'ZIZ  u.(x)  {i^(x')   =   e      "     p(x,x"') 


OCG  -1 ^r- '- 


where  H  =  p  /2ra  +  U  and  p(xjX')  is  the  Dirac  density  matrix. 

We  will  discuss  two  ways  for  evaluating  V(x',x)  and 
hence  M.(co').   Method  II  is  an  approximation  for  the  density 
operator  which  is  valid  when  the  potential  U(x)  is  slowly 
varying.   Method  I  couples  a  canonical  transformation  with 
a  variational  principle  for  a  Green's  function.   These  two 
methods  are  discussed  in  detail  in  the  next  two  chapters. 

In  order  to  solve  the  integral  equations  (2.55)  for  the 

12 
eigenvalues  we  can  make  use  of  a  variational  principle 

f  X.  (x)  V*(x,xSa3t)  4),  (x' )  dx  dx  • 

(5.6)   n.  (cD')=-^ ^ ^ 

^  J  ^i(x)  X^{t)    d^ 

in  which  M-.  (o)' )  is  stationary  with  respect  to  small  variations 
in  X(x)  and  *(x) . 

We  will  restrict  our  attention  to  the  lowest  eigenvalue 
which  we  will  denote  by  M-,  . 

In  order  to  choose  the  trial  functions  ^(x)  and  X(x) 
we  can  be  guided  by  the  following  considerations  I 

(1)  In  the  free  electron  gas  the  eigenfunctions  were  plane 
waves. 

(2)  The  lowest  eigenfunctions  should  be  spherically  symmetric. 


12 

Lovittj    W.    v..    Linear   Integral   Equations,    New  York, 

McGraw-Hill    {\^Wy. 
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(3)  Our  problem  Involves  a  .'   ,.,y  varying  external  i  .ial 

(4)  We  will  assume  the  external  potential  to  be  spherically 
symjnetric . 

Th  ■    r'e^  the  following  trial  functions  were  used  I 

X(x)  =  e-«l^l 

<l>(x')^  e-^l^'l 

where  a  and  p  are  complex  and  their  real  parts  are  non- 
negative  . 

Using  these  functions^  M-,  (co')  will  be  calculated 
explicitly  by  the  two  different  methods  referred  to  above. 
Then  from  M-(a)'  )  we  will  determine  whether  any  undamped 
collective  modes  exists  and  if  so,  at  what  frequency. 
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CHAPTER  VI 
The  Solutio^n  for  a  Dense  Electron  Gas  in _3n 
Ex t e rna  1  Po t ent_l_a^l^_-  Method  I 

1 .   Th 8  _S_olu t_i£n  _Expressed  in  Terms  of  a  Gree_n^'_s.  Function . 

Referring  back  to  Chapter  V^  we  see  that  to  evaluate 

M.(co' )  we  must  know  Y{x^  ,x) .      To  do  this  we  must  evaluate 
occ 

^__  u  (x)u.  (x«). 
i  occ 

We  can  rewrite  >   u. (x)u. (x')  in  the  following  form: 

i 

occ  •  J.TT 

YZ   u.  (x)ii^(x')   =   2<x  1  e"^^^  0(Ep-H)  [  x '> 
i 

where  t  -  t-t'^  and  H  is  the  hamiltonian  operator. 

The  factor  of  2  is  added  to  account  for  spin. 

Also, 

unocc  ,  i  -f-u         IV 

XI   u^(x')a^(x)   =   2<x>  1  e^^^  e(H-Ep)  ]  x> 
k 

As  stated  in  Chapter  V,  we  will  use  the  trial  functions 

<i)(x)    =   e-Pl^l 
Define  ^(rp)  to  tie  v  •  * 

q2 


v(r2,r)   = 

f{r^)      =  q^  J  -p r  d^^  . 

Ir  -r  1 
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After  substituting   into   equation    (5'6)    and   taking 
the   Fourier   transform  with   respect    to   t  we  have 


,3 


(6.1)      ^(-')    -^^-i^Jdt   e(-t)e--'t  Jd\dV 


it 


•[X(xi)^(x2)<xj   e^       P 


29{E^-R)\    x„> 


F 


it 


.    ^■ 


^      F         20(H-Ep)  I   x^>J 


,3 


+  %tPl^    r  dt    9(t)    e 
lii   O  IT     J  ^     ' 

it 


-icu't 


d^x-j_   d^g 


it 

h  ^^F 


(Ep-H) 


29(Ep-H)   1  x^> 


t  (E^-H) 
<x^l    e  ^      ^         2e(H-Ep)   I  X2>  } 


itEpA 


where  we   have  multiplied   and   divided   by   e 

The   first   term  is   the   direct   onej    the   second   is    the 

exchange   term. 

The   expression    M.(<^'  )   must  be   stationary  with   respect 
to  a  and   p. 

For  ease   in   calculation  ve   can   rewrite   equation    (6.1)    as 
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(6.2)      M.(aj'  )    =    Stat 


!«+§£     J  dt   e- 


ico't 


iiiSTT 


-i(E^-H)|  -i(H-E^)^ 

{9(t)    Tr    [X.    e  ^        ^    20(Ep-H)V'  e  ^    ^    29(H-Ep)] 


+  9(-t)    Tr    [X   e 


4-(H-Ep)  X-^H-E..) 

^  ^      20(E^-H)V  e^  ^    29(H-E^)]) 


Now, 


2<x^|    e 


^-(E,-H) 


(Ep-H)    I  X2> 


27ri 


^^iIz,+E^-hI^2^^ 


it 


zi 


dz- 


with  the  Z-,  integration  in  the  complex  z-,  plane  along  the 
following  path: 


z^+Ep-H 


N 


is  the  Green's  function  operator  for  a  system  whose 


Hamiltonian  is  H. 
it(H-E^) 


it. 


2<x^|e 


^'9(H-Eplx,>=  -^|iJ<X2l^-4-ZHl-l>e'  ' 


-^■ 


dz. 


The  path  of  integration  in  the  complex  Zp  plane  is  I 


'■>r^. 
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.o   equation    (6.2), 


(60)      M.(a)t)    .   Stat  iM§l^      r  dt   e"^^'^    f  dz,    dz^ 

it  it 

which  is    the    .  •  ■:      as 

,3 


(6.4)    M.(a3')    =   Stat   -^"-±^       J  dz      dz      Tr    [X  j-^-u   ^  yL    u] 


■  icu't      itz^A   -itz^A 


00 

•  [  J  e   —   -   e        ^  ^        dt 

0 

P    icu't      -itz^  A     itz    A 
+  J    e  e  e  dt]    . 

-00 

Combining  terms   and   integrating  with  respect   to   z   gives! 


(6.5)      l^(a)>)   =   Stat  ^°^-  r  dz      dz      {Tr    [X  t-J^-It   ^  ^^^-1?; 

.    r 1 +     1 1) 

"■Z^-Zp+O)'  Z,-Zp-CD'      J"' 

where 
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2 .   EvaluatlOQ  of  the  Green's  Function . 

For  an  arbitrary  potential^  we  do  not  know  the  explicit 
closed  form  for 


^^1  I  z,+E„-H  1  ^2^  > 


2m 


H  =  2_  +  V 


However^   we   can  make  use   of   a    variational   principle. 


—;jj     =      Stat    2G   -    G(z+W)G 
z+W  „  ^  ' 


so 


^'^ll^'+i"^l''2^     =    Stat  <x^l  (2G^-G-j_(z^+Ep-H)G^)|x2> 
IF  G-| 

The  problem  now    is    to  find   a    trial   G   . 

¥e   know    that   if  V  =   0,    then 

iik^lr^^-r^l 
<x    1g   |x   >     =      -  -^-2     --— — 

/r^l+Epj2H 

To   find  '\x    [  G-,  jxp)*     when  V"  ^  r     we   try   to  find   a 
canonical   transformation  from  the   system  where   V  7/  0  to 
a    system  where  V  =   0.      The    simplest   canonical   transformation 
to   try   is    one  which   transforms   only   the   spatial  part   and 
leaves    the  momentum  alone.      If  we   assume   V  to  be   spherically 
symmetric  we   can   restrict   the   transformation   to  f(|r[). 
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To  visualize  a  transformation  from  a  space  with  V(r)  to  one 
whare  V(r)  equals  zero,  we  can  imagine  a  new  space  where 
|r|  ■     retched  or  co:  whether  V  is 

repuls ;.  J  but  with  V(r)  •  ,    to  zero.   Also 

the  origin  and  the  point  at  infinity  must  he  unchanged  under 
f.   The  graphs  below  illustrate  these  general  features  that 
f(|r|)  must  possess. 


y 


V  repulsive 


y=r 


f(r)=y 


V  attractive 


Since J  in  this  new  space  V  =  0,  the  Green's  function 
is  given  by 

^±ikjf(r^)-f(r2)|   df    df 


<x^|G-|_|x2>  = 


m 


27rh 


dr,   dr, 


|f(r^)-f(r2)|   -n 

where  we  have  replaced  r  by  f (r)  in  the  G  for  V  =  0  and 
multiplied  by  the  Jacobian  of  the  transformation. 

Obviously,  there  are  many  functions  f(r)  with  these 
general  properties.   The  precise  f(r)  (which  may  still 
not  be  unique)  would  be  determined  from  the  specific  V 
in  question.   If  we  could  calculate  f(r)  directly  from  a 
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knowledge  of  the  potential  V(r)  the  Green's  function  would 
be  determined. 

However,  for  any  general  V  we  cannot  calculate  f (r) 
exactly.   It  is  not  necessary  to  do  this  however.   We  can 
choose  any  f(r)  with  the  general  properties  described  above 
and  insert  into  this  f(r)  parameters  to  be  determined  from 
the  variational  principle  for  G-,  .   From  now  on,  we  will 
understand  f  to  be  a  function  of  r  and  these  unknown 
parameters . 


where  the  unknown 


_  Q  I  -p 

For  example,  suppose  V  ■=  b  e 


We  can  use  f (r)  = 


iell' 


e^ 


±AA 


r  I 
parameter  e  is  complex  to  allow  for  a  phase  shift.   We  note 

that:    f(r)  ->  0  as  r  ->  0 

and     r(r)  — >  oo  as  r  — >  oo  . 

This  particular  form  for  f (r)  was  chosen  because  it 
satisfies  the  required  general  properties  and  somewhat 
resembles  the  form  of  V(r). 

Returning  to  the  variational  principle  for  G-,  we  shall 
use  as  a  trial  G-,  '. 

sin  (k,  If.  (r.  )-f.  (r^)! 
(6.7)   <xJgJx„>  ^  -  ^  A_lL^iAi_l_L2Z.L  .  G^Cx^^x^) 

(z,+E„)2m 


^-P 


F' 


i.2 
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We  have  c   l    .1  the  Jacobiai,  <         use  we  r.     ;  not  the 
exact  transformation  but  only  a  form  to  use  as  a  trial 
function.   The  sine  function  is  used  because  it  is  that 
linear  combination  of  e    /r  and  e    /r  which  is  finite  at 
the  origin. 

2 
=  J  d^x^  G-j_(x^,x^)(z^+Ep+  1^  V|-  V{r^)  G^^lx^.x^) 

We  want  to  evaluate 

<x^lG^(z^+Ep-H)G-j_lx2> 

2 
=  J  G^(x^,x^)(z+Ep+  1^  V^  _  v(r^))  G^lx^.x^)  d\^  . 

2 
2^  J  G(x^,x^)V  G(x-,jX2)d'^x^  can  be  simplified  by  using 

Green's  Theorem. 


J  div^  G(,^^)  grad^  G(3^2)  ^S  =  J  ^(l,3)S^^S^(5. 2)^^^ 

=  J  grad^  G^3_^^^.  grad^  G^^^^)  ^S  "^  /  ^(1.3)^  ^(3.2)^S 

Because  of  the  form  of  G,  the  surface  integral  (the 
surface  is  at  infinity)  goes  to  zero  and  we  are  left  with 

It  is  straightforward  to  evaluate  the  gradients  and 
to  take  the  dot  product.  We  leave  the  evaluation  of  the 
integral  to  numerical  methods. 
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t2 


Lit  US   denote  by   K-^(x-^,X2)    the   following  s\xm: 
(6.9)      K^{x^,x^)    =  J   G^(x-^,x^)     (z^   +   E^)G^(x^,x^r^   dr^d    cos    6^ 

-    2m  J  Vl^^l'^5^*Vl^^3^^2)    ^S    • 
Letting  R  =   f^(r^)    +  t^{rr,)    -   2f{T^)f{r^)    cos    {9^-9^), 


(6.10)  <x-]_l— :jr|-:H    '''2^   =    stat^(2G^(x-^,X2)-K-j_(x^,X2)) 
1      F 

sin(k  r1/^) 
=   Stat    (-  ^     f^ K   (x   ,x    )) 

Also, 

<^2U    +E   -h1^1>   =    '^'^f  <X2l2G2-G2(z2+V^^^2)l^l> 
2      F  2 


Try 

sin(kglf(rg)-f(r^)| 

<^2l^2l^l>    =   -2(^2^^l)    =    -   ;^  iTTr^)    -    f(r2)l 


with  ^__ 

y(z2+Ep)2m 


/ 


^2   -   v^  ^2 


As   before, 

(6.11)      <X2lG2(z2+Ep-H)G2lx^> 


2 
J  d^x      G2(x2,x    )    (z^+Ep  +  Ih  v2    -V(r^))G2(x^,x^) 


^3   ^2^^^2^^^5 
Define  K^{x^,Xj^)    as 
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^^         r.  ,  p 

(6.12)  Kglxg^x^)    =   -  ^  j   ^^G2(x2,x^).V^G2{x^,x-L)27rr^  dr^d  cos  0 
+  J  G^{x^,yi^){-    -t^'  ,)G^{x-,,x^)2ir  r'.,    dr     d   cos    9 

so  that 

(6.13)  <^2lF"-i:^"^l^^l^     =    stat^(2G2(x^,x^)-  K^{x^,x^)) 

2      F 

^       sin    (kgR^/^) 
=   stat^[    — -^     ^-5-^^ K2(x2,x^)] 

where   R  =-    .     ,rj_)    +  f^(r2)    -    2f(rj^)f(r2)    cos    {9^-9^)    . 

Expr    .  .    .10)  a^  :   '.   )  are  substituted  into 

/'  • 

Because  of  the  way  in  which  k,  and  k^  ^^^  defined, 

1/2      1/2 
t-.       rs  z,  '   and  z'      occur.   Therefore  we  must  define 

branch  cuts  in  the  z,  and  z^  planes  before  doing  the 

integrations  over  z,  and  Zp.   Mathematically,  we  are  not 

restricted  in  making  ths  branch  cut.        r  physical 

arguments  dictate  that  the  cut  correspond  to  the  bound 

states.   Actually,  if  we  make  the  cuts  in  any  other  way, 

we  find  that  ki.(cu'  )  is  identically  zero  which  is  not  a 

physical  solution. 

B'-^T '  w  -^re  diagrams  of  the  z,  and  Zp  planes  which 

indicate  the  paths  of  integration  and  the  branch  cuts. 
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— <— ^-, 

\ 

cut    ^^ 

) 

z 

cut 


We  now  refer  back  to  equation  (6.5) • 

-a   r- 


1 


r     e-Pl^i    .3 


jr-rpl 


d'^r    . 


3.     ^3, 


In   equation    (6.6)   we  have   to   integrate    over  d-^r.,    d-^rp  d'^r 
Let   d^rp   d''r  — >  hir  r   r^  drp  dr    [r-rpl    djr-rpj    d   cos    9 


(6.1^)    t^(a)')    =   Stat 


la+^ 


3 


apf        2Tr 


dz-^dZpr^dr^d   cos    i^^^^-^,^ ^,  ^-^,  ^^_ 


■,{r^,r   ,9      9    ) 


>-P   r 


r  dr  drpl  r-rpl  dj  r-rpl  d   '^°^ 


r-r^ 


where  g(r-,  j  rp,  0-,  ,  0p)  encompasses  all  the  terms  in  the 
integrand  which  are  functions  of  only  r,  and  r^',    i.e. 


(6.15)   g(r^,rp,9^,9p)  =  e'^'l  <x^\j-:^\x^X^^\~^^_-^\y.^> 


^1+^F- 


See  equation  (6.IO)  f  or  <(x,  ] — Tp  "-h"!  '^2^  ® 


nd 


2^rF 


see  equation  (6,1^)  for  ^^g^F +E^-h1'^1^  ' 
The  integration  over  [r-rpl,,  "^ >    ^^^  -"^p  ^^* 
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00  GO        r+Pp  CO  2        2 

Jdrgjdr        ;    ^Ir-rgl    +  Jdr^  J  dr     /     dlr-r^l 

0  0        r^-r 


0     rp 

r-r^ 

The  first 

term: 

oo 

^+^2 

;  rp  dr^ 

,^"  r  dr   f  ^ 

djr-r^l     |r-rpl    e  g(r-^. 


0  r^  r-r^ 

'-. 7-  d   cos    e 

Ir-r^l 

The  integration  over  dlr-r^]  gives  the  factor  2rp. 
We  now  have 

oo  00 

J   2r^  drg  d  cos  9^  s{rj^,r^,e^,9^)   J   e"^^  r  dr 

0  r^ 

oo 


-P  rg    rp   -p  r. 
2   ^  "*"  ^ 


P  -pr    ,       1    *^  "2    "2   '^  ^2 
I  e  *^   r  dr  =  -^  e       +  -^r-  e 


r  P 

^2 

The  second  term: 

o^  ^2  ^^2  /  ^  ^^  /  ^l^-^2l  l''-^2l  ^      g(ri.^2^^1^^2) 
0        0      r  -r 

d  cos  9      . 

r-rpi 

The   integration   over  d|r-rp|    equals    2r. 

oo  rp 

J  rp  drp  d   cos    9^   g{r-^,r^,9^,9^)   J  2r^  dr   e"^   ^ 

0  .  0  ^ 

To  "*"    "'^2      2  "^   ^2  -Br 

r  dr   2r2   e-^-     =      -  !l -^-  -  --^-l- ie        ^  _^  ^ 
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Defining; 

-Pr2      2        ,          -P^2 
2e             Pp        4r^e 
gplr^.r^^QjL^S^    =    g{r^,r^,9^,9^)    ( ^ ^^ 


(6.15)    ^x(a..)   =   stat^       i«gi-     J  dz^  /  dz^  /  dr^   r^ 

^  0 

IT  TT 


•   J    Sin   9^de^  J    sin   92^^2 

-TT  -IT 

GO 

•[   J    2r2g^(r^,r2,0^,92^^1^^2^^^    ^^2 
0 

00 

+     J  ^2   ^^2    S2^^1^^2^^1^^2^^1^^2^^^ ]    • 
0 

We   leave   the   remaining   integrations    to  be   done 

numerically.      Then  we  must   find   the   stationary   value  ^by 

varying  a,    ^,    and   the   parameters    in   the   transformation  f. 
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CHAPTER  VII 
The  Solution  for  a  Dense  Electron  Gas 
in  an  External  Potential  -  Method  II 

1 .   An_  Approximation  for  the  Modified  Interaction  for 
:         Varying  External  Potential . 

We  saw  in  Chapter  V  that  the  solution  of  our  problem 
depends  upon  calculating  the  eigenvalue  M-  (co'  )  which  was 
approximated  by  equation  (5.6): 

J  X(x)  V*(x,x'^a)')  <})(x')  dx  dx"' 


(5.6)   M.-,(a)')  = 


<D(x)  x(x)  dx 


where^    for  a    spherically   symjnetric   potential  we  decided 
to  use 


X(x)   =   e-«l^l 


<t)(x'  )  =   e 


-pIx'I 


,* 


and  V  is  the  modified  interaction.   This  expression  for 
M--,  (cD'  )  niust  be  stationary  with  respect  to  the  complex 
parameters  a  and  p.   The  modified  interaction  V*  is  given  by 

(7.1)   V*(xSx)  =   J  V(xSx")  -^if^^^     dx" 

|x"-x  I 

It  was  shown  in  Chapter  V  that^  in  the  random  phase 
approximation J  V  is  given  by 
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UnOCC  OGG 


(5.4)   V(xSx)  =  ^   [{   Y7-  '\i^')    \i^^)   ^'\i^^)    u.(x-))9(t'-t) 

K  1 

UnOGC         _         OCG 

+  0(t-t')  ^  u  (x)  u.(x-)  ^_   u  (x-)  \(x))] 
i  k 

where  ^(x)  =  1  ^   x  >  0 

0  ,   X  <  0 
^  ,      X  =   0 

and  Y. '"-i-(^)  ^.(x')  :=  e~   ^  ~  '  ^    p(xjX') 

2 
where  H  =  |—  +  U 

and  p(xjX')  is  the  Dirac  density  matrix. 

If  the  potential  U(x)  is  slowly  varying,  we  can 
approximate  p(x^x')  by  I 


00  +i£ 

(7.2)   p(x,x')  =  -  -S^-      r    do)  f  d^p 

h-^27ri  ^  ,  .  ^   "^ 

-  OO  +1  £ 


i?.(5?-^")A   -i^[Ep-  1^  +U(x  )]+  ^^4l^U(x)|2 


CO 


,  ^      X  +  x' 

where         x  =  — ^^ — 


This  approximation  is  a  corrected  Thomas-Fermi  model. 
If  Vu(x)  is  zero,  then  (7.2)  is  just  the  Thomas-Fermi 
result.   A  complete  discussion  of  how  this  approximation 
is  arrived  at  is  to  he  found  in  Appendix  A. 
To  find  e-i^(*-^')/^p(^/x^'): 
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(7.5) 


p(x,x')    =    9(E„-H)    5(x-x') 


oo+i£     -ia)(E„-H) 


-^  r 

27ri     J 
-oo+ie 


cu 


■■-?. ) 


so  that 


-iH{t-t')A     _  _^  -iH(t-t')A      n 

(7.^)      e  p(x,x«)   =   -e  ^ 


oo+ie     -ia)(Ep-H) 

e  —  (jcu  5(x-x'  ). 


CD 


-oo+ie 


Let  u  =   00  -    (t-f  )  A  and      co  -  u   +    (t-t '  )  A.      Then 

-iH(t-t')A  ,    ,  -IE    (t-t')A 

(7.5)      e  p(x,x')    =   e 

CD+ie   -iu(Ep-H) 
p        e 
J         U4 
D+ie 

For  slowly     varying  u(x)   we   can  approximate    (7.5)   t)y 


-oo+ie 


-iH(t-t') A 
(7.6)    e  P(x,x') 


2  -iEp(t-t')A 
e 


h^2iri 


QD+ie                          i  g.(x-x')      -ia3(E^-  ^  -U(x) 
/     d^o  J  d^p  [  ^ 


2 
)_ 

'F"   2m 


-oo+ie 


OD+(t-t'  )  A 

ia>5ii^|vu(x)|^   /24m 


] 


Also, 

unocc  -iH(t'-t)A  _,  ^     ^ 

(7.7)  y_:  Uj^(x')    aj^(x)    =   e  [5(x'-x)    -   p(xSx)] 
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and   if  U(x)    is    slowly  varying 

-iH(t'-t)/ti      ^     .^  +iH(t-t')A        -^     _> 

(7.8)    e  5{x'-x)   -   e  6(x'-x) 

i  ^~    t&m-  +U(x)]+ii:^;l^-^-lvu(x)l2      i?.(^.-^)A 


h^ 


For  the  time  being  consider  only  the  first  term  of  V(x'^x) 
since  the  second  term  differs  only  by  an  interchange  of 
X  and  x' . 

_  unocc       _     occ 

(7.9)  V(xSx)  =  ^  [e(f-t)  YZ   Ui,(x')  ^^^(x)  YZ   u.  (x)u.  (x-)] 

_  -  — (t-t  '  ) 

(7.10)  V(x',x)  =  3^  e(t--t)  e   ^       p(x,x') 

+1  H[t-t_:l  ^.  H(t-t') 

•  [e      ^    6(x'-x)  -  e      ^    p(xSx)] 

Writing 

-oo+ie 

using  the  aforementioned  approximations  (7.6)  and  (7.8) , 
and  taking  the  Fourier  transform  of  V  i.e. 

i.Je        ^    ,^(t_^)  V(x«,x)  , 

and  also  letting  x  =  r.^   x'  =  r   we  get  the  following 
expression  for  V{r^,r^  ,od^  )l 
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2-1-    ^    -    .^5-^2   J 


(7.11)    V{r^,r,,G,')    = 

-^^(V   2H  -U( 


•1 


R 


•i^.(^) 


''     -'^p).'-^\^.i'-^)\^ 


/[^l  +  (^)]a3 


i(t-t.)[?|  .u(^).  lit^'lM|vu(^'-^)l2 


i^(^)       -103- M^) 


t-t-x   ^5.     ^3, 


icu^    dCD  d(^^)    d^p-^    d^Pg 


iP]_-(r-j_-r2)      ip^- {r^-v-^) 


co'+ie       oo 


"^  -oo+ie      -oo 

ip-^-lr-^-r^)        ip2-(r2-r^) 


[e 


-io)  cct   +ioa'tiOD. 


^JXoj^+oj^)  ] 


-103. 


1 L  F   2m 


p^    ^1+^P    i^^^,    ^1+^p  ,o 


^1 
[Ep-  —  -U( 


m 


00    00 

-1    J,^^  J  ^^Pi  'i^P2  n  ^^2 

in  TT      CD»  -00 


r   iPi'(ri-^2   iP2'^^2"^l) 


-^•[e 


[e 


P?    r  +r     io^^ix^     r, +r„   „ 
-ntV  2E  -U(-%^)]+  ^^-  |VU(-V^)|2 


•  e 

We  assume  that  v( | r, -r^\ ) ,    the  interaction  between  the 

^        2 

particles,  is  the  coulomb  interaction,  v(lr-,-rpj)  =  -i — ^ . 

I^l"^2l 
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(7.12)   v'^Cr  ,r,a)')=  [—3- -v(r   r^,cu')  dr 

I  2   I 

and 

-ar   -Pr   ^ 
e     e    V  (r-,,r,cD'}  d  r^  d'^r 
(7.1J)    t.,{co-)==^ r:-a^-Pr  .3. 


/ 


e    ^     d  r 


2 .   A  Sample  Term  of  the  Solution  Expression  as  a  Many 

Dimensional  Integral . 

P  ft        ^  ^^1  (^'  )    Stt 

Call   h  (^'  )        e"^""   e"   ^"^  d^r  =  M-,  (o)'  )    =  -± ^ 

_  (ct+P)^ 

¥e   now  lA/ant   to   calculate    |i.-j^  ( cd '  ) . 

Since  the  different  terms  in  V  all  involve  the  same 

integrations  in  r^r.,  ^r^j  p,  and  Pp   just  look  at  the  last 

term  to  do  these  integrals.   Let 

-ar      p  ^    <l)(k  )   ik.-r 

(2Tr)- 


(7.14)   e   ^  =  J  d\  ^^-  e 


-pr     r,5;.  ^^^2)   i^2-^ 

*(^i)  =  T-f-^  ^    ^^^2)  =  T-^\2  • 


Also,  make  the  following  sutstitutions : 


(7.15) 


^1+^2 
2 

=   Pi    > 

Pi   +  9% 

— » 
=      ^1 

^1   ^2 
2 

-h  > 

Pi   -   P^ 

-      ^2 

— >        — > 
rp-   r 

=      92' 

— >         — ) 
Pi   -   P5 

-    P  2   = 

The  Jacohian  of  this  transformation  is  2. 
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The   las 


)    is 


oo 


)       M-il^')    - 


oo 

r 

•   J 

CO' 


dco 


-00 


ih    TT 


ioxD  li-icD'hO) 


,  ,,^  d^k^d^k2 


*{l^l)    ^(kg) 


CD 


(^1+^2)  IP2I       (2^)^(2^)^ 


ik-]_(p-L+p^)      ik^Cp-^-p-^-pg)      ipj^-2p    A      -iP2-2p    A 

•        p  -^p  -^p  -^p  -^ 


•    e 


-i^2tV  2H-u(Pi)]+^fijrl^u(p,)l2] 


The  integral  over  pp  is 
(7.17) 


J  TT  ^2  ^  ;;2 


the 


The  integrals  over  p,  ;,  p^  and  p^  combine  to 

following 

^k^.,    ,3/2  ,2      it^lk^-k2lS^2  /8n^(^+^2) 
(7.18)    ^^if )    \yg   e 

and  the  last  term  of  jl-,  (cd'  )  becomes 

°°f  ^   ?^  d^^  da3_  d^p,  a\    d\. 

(7.19)    J  do.  r  ji_^ — ^:i.i — 2 

co'+ie  -00      ^  12' 


^   icugh-io)- 1.0)2  ^  ^_8(2m)_5_{V 


i^'^h^;^"    ^  (k2l2)2  (^^-2;^ 
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il?^-p^   ii?2*Pi    i|k^-k2l^^^^l'-^2  ^^'^^^+"^2^ 


4 


e 


-icD2(Ep-U(p^))+  -^^  lVU(p^)r 
e 


3 .   A  Variational  Principle  Is  Introduced  to  Simplify 
the  Integration. 

We  now  consider  the  integration  over   p  .   If  the  form 
of  the  potential  U(p-,)  is  particularly  simple  (for  example^ 
linear  in  p-,  )  then  the  integration  over  p,  is  straightforward 
and  not  difficult.   However,  for  many  more  complicated  forms 
of  U(p-,  )  the  Integration  can  no  longer  be  performed  in 
closed  form.   Also,  it  is  preferable  to  continue  the  problem 
without  having  to  make  a  specific  choice  for  U(p-,  )   at  this 
point  in  the  calculation  so  that  the  final  solution  would 
be  valid  for  any  general  U(p-,')  consistent  with  the  initial 
assumptions  of  spherical  symmetry  and  slow  variation  of  U 
with  p-]_.  To  accomplish  this,  use  was  made  of  the  following 
variational  principle! 

(7.20)   J  e"^^^^  W(x)  dx  -  J  f  (x)a)(x)  dx  -  J  u(x)a)(x)f  (x)  dx 


f(x)  In  f(x)  a)(x)  dx  . 
where  the  expression  is  stationary  against  variations  in  f, 
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'i-  -.-jw  the  validity  of  this  varigtional  principle, 

take  the  variation  of  the  right  hand  side  o-  ;   ;t  it  equax 
to  zero,  i.e. 

0  =   [J  f(x)cD(x)  dx  -  J  u(x)co(x)f  (x)  dx  -  J  f  (x)  In  f  (x)a3(x)dx 


0  =  J  5f  (x)co(x)  dx  -  J  u(x)cd(x)  5  f  (x)  dx  - 

'5  f(x)  In  f(x)  cd(x)  dx  -  J5  f (x)a)(x)  dx 


u(x)'-d(x)  5  f  (x)  dx  -  f  5  f  (x)  In  f  (x)  a3(x) 


dx 


0  =  -  J  o^{x)  B  f  (x)  [u(x)  +  In  f  (x)  ]  dx 


In  f(x)  =  -u(x), 

f(x)  =  e~  ^  '   at  the  stationary  point  of  the  right 
hand  side.   With  this  value  for  f(x)  the  right  hand  side  is 
J  e-^(^)'a3(x)  dx  -  J  u(x )a3(x )e-^(^)  dx  +  J  e-''(^)u(x)'i3(x)  dx 


=  J 


,-u(x) 


a)(x)  dx  which  is  the  same  as  the  left  hand  side. 


r  -^(pi)  ,  ^  3 

We  want  to  ca^lculate^   e      "^(P-i  }  d  ( 
i(i?,+ic  ).p' 
where  a)(  p  )  =  e 

i(a^-K.i|)t2lvu(pj|2 
u  -  -  i(co^+co2)U(p^)  -    -"   ^ 


24 


m 


Assume  f  =  f  ( [  p^| ,  p.^  .  ..p^) 

where  the  p-,...p.   are  some  set  of  i  parameters  to  be  varied, 

The  spherical  symmetry  of  the  potential  has  been  made  use  of 
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by  assaming  f  depends  on  j  p.  j  and  not  on  p  . 

Consider 

(7.21)   -  J  f(lp^U  P-^...  P^)  m  f(lp-i_U  P^  ...p^)a)(p^)d^p^=I. 

Ig  =  -  J  f(lp-Lh  Pi  ..-Pi)  In  f(lp^Up^...p^) 

il^l  -Ip Jcos  9 
•  e  2Tr  d  cos  9    [  p,  |   dp^  , 

— >   — >    — ¥ 

where   k  --^  k-,  +  kp  . 


I2  =  -  2Tr 


f(lPlh  Pi.-.P^)  In  f(lp-j_l,p^...p^)|p^l  dp^ 

ilkllp  1     -i|k|lp  1 
.(e      ^   -  e       ^  )  /ik|p3_l  . 


In  a  free  gas  case^  the  plasma  mode  occurs  in  the 
region  [kj  <<  ]kjp].   To  simplify  the  calculation  (especially 
the  subsequent  integration   over  k  +  kp)   we  shall 
make  the  same  assumption  hare.   Since  the  purpose  of  the 
calculation  is  to  study  the  plasma  region  and  since  a  plasma 
mode  is  basically  a  long  range  correlation,  this  should  be 
a  valid  assumption.   If  our  purpose  were  to  study  the 
correlation  energy,  short  range  effects  are  important  and 
the  assumption  jk]  <<  lk„)  would  no  longer  be  valid. 
Furthermore,  if  we  are  interested  in  a  system  such  as  an 
atom  1  p.,  ]  is  never  much  larger  than  10~   cm.   Therefore, 
we  can  expand  the  exponential  term  and  keep  only  the  first 
nonvanishing  term; 


i   k|    P-,           -i   k 
e                   -   e 

Pl 

=    2    . 

i   k     p^ 
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A-tually,    we  will  presently   show    that   when   the   specific 
potential    is    introduced   and   all   the   integrations    over   p, 
performed   before  making  any   approximation,    we   can  get   th:? 
same    result   without  m:J..^..^^    ...y   assumptions   about    p,  . 

With   the   above  approximation     we  have 

{7.22)      ig  =   -  4Tr  J  fdp-^l,!  |.  ..P^)  In  f  ( ]  p^l  ^p-^.  .  .p.  )p^  dp^, 

Consider  -        ^(pi  )<^(p-i  )f  (p-,  )    d'^p-,.  The    :  term  is 

(7.23)      I3   =   2Tn{ax^+^n^)    "  •  (!p^1)    e  ^f(lpil)   d   co.         ,^:, 

and  as   kA-o   ^^  ^, 

.       )      1^  =   +  Hiri^o,^+<o^)   J  u(1pJ)  f(lp-Ll)    pI  dp3_    . 
Let 


)      14  = 


^      h- 


27ri(a^+a^)ii2      ,    ,         ,       ,    ,p  , 
4^m— -     i    I^U(lPil)l      ^(Pi)    f(Pi)    d\ 

1^^—     J    |vu(1pJ)|2   f(lpj)    p^  dp^    . 


(7.27)  J  fdpil  '"(Pi))  d^p^  =  i^     ; 

(7.28)  I^  =  J  fdpje""^'^!)   d\    . 


:  potential  udp^^l)  ro,   f(lp^|) 

;    ^^    .....    oecause   at   the   stationary  point   f(x)    is   given 
.,,rv) 
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Wlien  f(lp-,l)  — >  1,  I-,  appears  to  be  singular  when 
k  — >  0.   Actually  there  is  no  singularity  after  the 
integration  over  k  is  completed.   However^  it  is  convenient 
to  subtract  out  this  point  and  treat  it  separately   i.e. 


(7.29)   I^  =  J  (fdpil  -1)  e    1  d^p-|_  +  J  e    ^  d^ 


Call  only  the  first  term  I, .   Call  the  second  term  I  . 

(7.30)  I^  =  {2-irp   5(k) 

(7.31)  ^1  =  ^^  J  (^(Ipil)  -l)Pi  ^Pi  • 

Now  suppose  the  specific  potential  we  are  interested  in 
has  the  form  U(|p-||)=  b  e      .   We  now  need  a  trial  f(|p-,l). 

Referring  back  to  the  variational  principle  we  can  see  that 

-u(|pil) 
the  correct  f  would  be  e        .We  would  like  the  trial  f 

-u(pn  ) 
to  be  as  close  to  e       as  possible  but  easy  enough  to 

handle  mathematically.   A  good  choice  for  f,  with  the  specific 

potential  U(p. )  mentioned  above  is 

where  A  and  Y  ai^s  the  parameters  to  be  varied. 

I2  =  -  47r  /  f(lPil)  in  f(|pil)Pi  dpj_  , 
-Ae-'^lpl 


w 


ith  f(lpj)   =   e 

-tIpiI  _y1p  1  -'^IpiI-aId  1 

I        r    -A  '^^                'iPii  ?              ,        p  _A^        -^     '^l  Pi  1  p 

1^=  hiTA  j    e  ^   ^               e  Pi   dp^  =  hwA  Je   ^^             e             p^dp^ 


A=Y 
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0 


Let  X  =  Yp-| ,   dx  =  dp,  Y. 

2    -  A 
(7.52)     I   =   ^  [  -^  [A   '^  T(A  ,A)]] 

'     OA  ^  A=T 

where  y(—  ^A)  is  the  incomplete  gamma  function.   In  the 
same  way  we  get 

4vri(co  +a)^)b   p^2   -  ^   . 

(7.33)  1.   = ^7-^-  [^  [A  ^Y(^,A)]] 

■^        "^       dA^         ^  A=a 

(7.34)  I^  ^  1-^ [-'^-^  [A   ^Y(^.A)]] 

24  m  Y       ^A""         "^      A=2a 

2    -  A 

(7.35)  li  -  T^  [A  [A  "^  G(A  ^A)]] 

A 
where  G(^  ,A)  =  J  x^'^^  ^"^  dx 

0        "    ,   1 
Using  the  potential  U(lp|)  =  b  e      ,  the  p, 

integration  will  now  be  done  before  any  assumptions  about 

k  are  made  in  order  to  show  that  no  assumptions  about  the 

size  of  p,  are  necessary. 

We  want  to  evaluate 

r  "^(Pi)       ^ 
J  e       (o(p-^)  d^p^ 


i(k,+kp)-p-, 
where  cD(p  )  =  e 
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i(a^+a|)t2 


U 


(Pl) 


=     b   e 


-ap- 


VU(p-^)l^        =     a%^  e  ^    . 

Say        =   e  .      Than  A   and   y  are  parameters   which 

— >       — »  — > 

will  be   varied.      Let   k  =   k-,    +  kp. 


-   /  u(p^)a3(p^)f(p^)    d^p^   = 


2F(co^+a)2)b    r^        -A/T  1 


T|kl 


^(A  7(^   ,A) 


J>f=a-i|k] 


+ 


Tiki        1  |A=a+i|k| 


^^""^^'^'''^'^'^-[.-/-n^A)] 


24  mYlk 
"•"  2^  mYl    k 


^-^   [A~"''    y{^  ,A)] 


A=2a-i   k 


A-=2a+i   k 


J  fdP;^!)   m  f(p^)    e  1   d^p^ 


^  fR{^t'^"''^^(^^)^] 


A^Y-ijkl 


_   2T:lAr^        -A/T  ,  ,A      ,s.\ 


Y|k 
f(Pl)    e 


5a 


[A-"/''   T(^   .A)] 


A=Y+i|k| 


^   d^p^      =     J    (f(p^)-l)    e  ^  d^Pi+Je  ^d^p 


1' 


The  first  term  is 
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A 


A/T 


where       ^  x  (e   ^^-1)   ^     is      G(^   ,A). 

tj  A.  [ 

0 
The   second   term  is   just    (2F)'^5(k)    =   Ij-. 

In   order   to  be   able   to  perform  the   integrations    over 

k-j^   and   k^   it  was   necessary  to  make   the   following  approximation, 

Let 


1^    [A-^^/^  G(^  ,A)]      _       =   g(k)    . 


A=+ik 


Let 


A'    =    ^ 


^F    ' 


y '    =     ^ 


k^'' 


-A'/Y' 


"^v'^^K. 


s(ir-) 


Assume    \r—\    <<   1, 
^F 


^F 


Call  kAir  =   k'j 


g(0)    +— ^ 


dg 


<i(^) 


_k 

k 


^F  ^0 


F 


dg       _     ^g       dA' 


dk" 


^i^rx 


Sa        dk'  Sa' 


(+i)    . 


The   terms  with  g(0)   all   cancel. 


-75- 


>2 


A' 


dA'"^  "^  A'=0    F 


I.    - 


ik-p.       ^ 
-1   -  J     (f  (Pi)-l)    ^  <^   P-^    "- 

Also, 

,5 


A' 


^        Yk^     5a-2 


[A      '^'    G(^   ,A)]] 


A'=0 


-J    u(p^)a3(p^)f(p^)    d^p-^ 


AJ_ 

T'    ../A' 


Y    i^T 


Sa'^  "^  jA'=aA, 


+ 


47ri(a^+a|)h^a^b^    r  ^2 


24m  Y  Kp  l^A' 


[A      ^'    T(^    ,A)]' 

■^  -*  A  '  =    2a  Ax 


15   +   14- 


ik'p. 
f(p-j_)    m  f(p^)    e 


47r  A 
,  2 


r  a= 


Sa' 

A 


"^  Ja'=yAt 


=  I. 


Since   A'    =  i^^  ,    Y'    =  — 


^F    ' 


4  A  [    a 

^    Ua 


2         -A/Y 
Q     [A 


Y(f   ,A)]' 


A=Y 


The  same  is  true  for  I,,  I  and  Ik.   The  agreement 
with  the  previously  derived  I-,,  !„,  1^,    Ik  and  1^   is  obvious. 
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4 .   The  Sample  Term  for  the  Solution  _i_^  Reduced  _t_o_ _a_ 
Two  Dimensional  Integral. 

Consider  equation  (,._.j;.   The  integration  over  p, 
was  discussed  in  the  last  sect"    ^nd  the  results  ?j  ■ 
in  the  factors  I-  1^,    1    ,    1^,    and  I  .   Suppose  we  consider 
the  in.,„3ijtion  over  k,  and   Ic„.  for  the  terms  with  I,,  I  , 
1     i  Ih,  but  not  I_.   By  using  contour  integration  and 
dec  ,    Ltion  into  partial  fractions,  these  integrals  were 
carried  out  in  closed  form.   Although  quite  lengthy,  this 

■  ion  was  str'  '    forward  and  the  details  will  be 
:    .  in  Appendix  B. 

If  we  substitute  the  factors  (7. 20),  B„   (see  Appendix  B), 
(  .  "] ,    (7.24),  (7-26)  and  (70l)  into  the  expression  (7.19) 

we  neglect  for  the  moment  the  contribution  from  the 
t-;rm  with  I^-,  then  the  last  term  of  P.(cjd' )  is: 

00    00    00 


(I^+I^+I^+I^)  e 


00(0)^+0)2)-^/ 


2       ..2  2 

(1  -  it'^''/^)) 


^m 


2C^iti^xTr(a-(^)   -ixii^(a+p)'^/8m  ,-^ 

--T„ <=  (l-ilMa-HB)  /^  ))  - 
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2iC^(a+p)hV    ^        ^ 


52m 


2 


-IX 


2/  ,^2 


7r(a+^)    e-i-^'(«+f^)'/8-    (l    -  ^(Ma+p) /^   ))] 


where   x 


^^2 

CD-  +C0^ 


A, 


The   integration   over  k-,    and   k„   for  the   term  vjith   I_ 

must  now   be   done. 

Remembering  that    I^.   =    (2Tr)^6  (k-, +lCp)      we  must   evaluate: 

i|k    -k    lt2a.^C02/^m(a3  H-ri)    ) 
^  ^  P      ^        6(k  +k    )a   p    e 

(7.57)    (27r)5Jd\     J    d\,  —  ^     ^g  g^g    ^^_,    ^    ^g^g^_g- 


-oo 


=   a  p(27r)^ 


-00 

„     d'^k-,    e 


(kf   +  a^)^    (1,2    +  p2)2j^. 


oo    .3,.      _ilkil^^S^2/2^(^+^2) 


^1 


-CB 

oo 


"(k2-K.2)2    (k2+p2)^2 


2.  2 


=   32Tr"^a  p  J     d|k^l 
0 


1^1^ f    (k2+p2)2    - 


Making  use  again  of  the  method  of  partial  fractions! 


E. 


E. 


(k24<x2)2(k2+p2^2  -  (1,2^2^2 


+ 


(k^+a^) 


+ 


+ 


(k2+p2)2    (^2^p2^ 


E 


E,  = 


-2 


-2 


2    /„2  .2^2  '   "1  ""  /^2   2^  ^   "'2 


(P  -a")' 


O  -«  )- 


/  2  „2x2  ^  "1"  ,  2  „20 
(a  -p  )         (a  -p  )^ 


The  integration  was  then  completed  by  using  standard 
integral  tables. 
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We  substitute  the  integrated ■.;.  .into  the  expression 

(7.19)  ■     t  the  correction  to  p.(''-D' )  due  to  the  term  with 

I^.   This  correction  is: 
0 


(7.38)  ---^4^4,^  r ,,  J  ,^^  J  ,„^ 


00        00  00 

dcu. 


CD'      -00        -00 


iCD^cjoh-io)' hoo2-iEp(ax^-HJ02) 

1 


1  iii^x 

/q2      2>2        2m 


11   -  I  (a^  /4l  )  1  Id  « 


-ia^ii^x/2m 


1  ih   X      v^  y/^ 

(0^-a^)^     a^2m     2!i/:Tr 


+ 


v-2^ti  -  i(ak  y^ )  ]  ^ 


-ia^ii^x/2in 


a    (6   -a    ) 

__1 

/    2  «2x2        2m 
(a   -p    ) 


in    X 


-IX     \  1      TT 


[1  -i(ei^  ^  )] 


-ip^ii^x/2m 


2m   ' J    2p 


+ 


ih  x        /f     ■,/^_ 

(a2-p2)2   ^2    2jn   2t.      ^^^^ 


-ip\^x/2m 


-iii^a^x/2m 


(aV)'   2P 


Ib  ti-i(B  W^a)]  - 


-iii^p^x/2m 
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5.   The  Complete  General  Solution  Expressed  as  a 
Two  Dimensional  Integral. 

If  we  refer  back  to  equations  (7.9)  and  (7. 16);,  we 
see  that  we  have  calculated  only  part  of  the  contribution 
to  jI(cD').   However,  as  was  stated  before,  the  additional 
terms  are  so  similar  to  the  ones  we  have  calculated  that 
we  can  evaluate  their  contributions  without  redoing  the 
calculation  for  each  term.   Consider  first  all  the 
contributions  to  the  direct  term  for  V.   We  notice  that 

the  integrations  in  r,    ^-i^^p'    ^1  ^^'^  ^2  ^^^   ^^^  ^^^   same. 

The  contribution  to  iI(co' )  that  was  calculated 
corresponds  to  the  last  term  in  V.   This  term  contains  the 
following  integral  with  respect  to  co: 

00       iCD  tO) 

f  ~^^^~^aj ^^'^   ^^^^    equals   -  E^(  ito^  ii  oj' ) 

03' 

where  E.  ia  the  exponential  integral. 

Consider  the  second  term  in  V{r^,    r-,,co').   This  term 

CO'     -icQ^coh 

e  dco  /  03  which   equals 


contains   the   factor: 

-00 


,  -00 

^        -10)^0^1  „  -ico-,a)  ii 


■f   (e  d'i)  /o))-   iTT        ^{^)    dco  e 

-00 

where  i  stands  for  the  principal  part,  and 
J-  e   -^   /cD  ^   -  E^(-ia)' ii  01^)  . 


-00 

Consider  first  the  contribution  from  the  principal  part, 
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If  we  r.^^i.^e   od'  by  -^JO' ,  we  see  that  this  term  looks 
exactly  like  t.'.    :'      '  last  term  exc       ■    Inus  sign 
and  the  fact  t^     .        of  ox^  and  oo  are  interchanged. 
Since  the  paths  of  integration  for  cu  and  co  are  the  same, 
we  can  u.'      ■    lange  their  names  without  changing  the     j 
of  the  XI.  ,    nee,  this  term  is  the  same  t 

term  with  the  following  changes!  (l)  multiply  by  -1, 
(2)  let  CO'  ->  -a^<  . 

By  the  same  reasoning,  the  term  with  6(^)  is  the  same 
as  the  last  term  with  the  following  changes  I  (l)  multiply 
by  l:r.    ),    let  co'  — >  -o)' ,  (3)  omit  the  exponential  integral 
factor. 

In  the  first  term  of  V  ,(t-t')A  plays  the  role  of  co 
;et  the  same  contribution  to  p.,  (00'  )  as  from  the  last 
term  with  -  E.[itia)  od']  replaced  by  2Tr  i9(-biCD  ). 

Now  consider  the  exchange  term,  the  second  term  in  (5.^). 
Again  the  calculations  are  very  similar  and  the  results  can 
be  obtained  by  simple  changes  in  the  direct  term. 

In  the  first  term  of  V  exchange,   (t'-t)A  plays  the 
role  of  00    This  term  is  exactly  the  same  as  the  first  term 
of  V  direct  with  'X)'  — >-<'X)'. 

The  results  for  the  other  terms  are  such  that  all  the 
terms  involving  the  exponential  integral  fimction  cancel 
those  of  the  direct  term  and  the  part  arising  from  the  point 
contribution  -i5(a))  is  exactly  the  same. 
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Let  1{(X)^  ,ad  ,CD  )    be  the  sum  of  the  integrands  of  the 
expressions  (7.36)  and  (7.38).   Denote  the  sum  of  the 
contributions  to  p.(co' )  from  all  the  terms  by  i^(^' )t;otal* 
Then 
(7.39)  M.(a.')total  =  1^1  ^^2  ^^^^  K-^'^^  ^g) 

Note  that  the  integration  over  dco  has  been  accounted 
for  in  this  expression.   From  no^^J  on  the  subscript  total 
will  be  omitted  and  il(a)'  )  will  be  understood  to  be  the 
total  value. 

6 .   This  Method  Applied  to  a  Uniform  Electron  Gas  -  The 
Solution  for  a  Uniform  Electron  Gas   Expressed  as  a 
Two  Dimensional  Integral. 
It  is  worthwhile  at  this  point  to  repeat  the  calculation 
of  P-(co'  )  for  the  free  electron  gas  case  using  the  same 
techniques  and  orders  of  operations  as  for  the  inhomogeneous 
case.   Of  course  J  this  is  a  roundabout  and  clumsy  way  to 
handle  the  free  electron  gas  case.   However,  there  are  two 
important  purposes  in  doing  this.   First,  it  is  a  convenient 
way  of  checking  p.(to',cLL  co  )  for  the  inhomogeneous  case  since 
this  must  reduce  to  il(co',ai^  cu  )  as  the  potential  goes  to  zero, 
or  to  a  constant.   Second,  since  the  integrations  in  oi^  and 
oi  are  difficult,  we  can  be  guided  in  the  method  by 
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performing  them  in  the  free  gas  case  where  we  already  know 

the  answer. 

Refer  now  to  the  e:-;.       ■•    for  .  .  j  for  i^(od' ) . 

In  the  free  gas  case  wa  .....'^...iy  know  the  values   for 

X(r)  and  0(r).  _^  _^ 

ik*r' 
X(r  )  =   e 

-^  — > 
-ik«r ' 

^('r')   =   e 

Also,  the  expression  (7.2)  for  p(xjx:')  with  U(x)  and 
VU(x)  both  equal  to  zero  is   the  exact  density  operator 
for  the  free  gas. 

We  repeat  all  the  steps  as  before. 

The  expression  (7. 16)  is  replaced  by 

00  00 

dco 


(7.^0)    |x(fjo' )    =   |a(a3' )    =  J  d^pj_  d^pg  do^  dco^  d^p^  <3?p^  ^^pj' 
-2q' 


-00  ^' 

,2         /--2^-i^-'^-2      _^_     il?.(?^-f?^)    -ii^-(p,-p3-p2) 

ip    •2p    A      -iP2*2p3A      -ifi3   (E   -p^/2m)-icu  (E  -p2/2m) 
•    e  e  e  ] 

Notice  that  this  can  be  obtained  from  the  general  expression 
simply  by  the  following  substitutions  I 

<D(k^)     V'lkg) 


27r^      ( 27r )  • 


=  1 


k  —  k-|  —  —  kp 

u(p-l)  =  0 
vu(p^)  =  0 
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Continuing  the  calculation  we  find  for  M-   ^(^')^ 

ik  ix  CO-  CO 


F 
1^^2' 


i^/"^h^  /f  k^     J     1    ^  (m  +CO^p/^ 


ico'iicb2  -ito'iKXi 

[2fri    e  -2iri   e  9(-iia)   ) 


icD'iico 
27ri   e  ^   e  ( -   0^2 )  ] 


7 .   The  Limit  of  the  General  Solution  when  the  External 
Potential  Goes  to  Zero. 

We  now  want  to  show  that  the  general  expression  (7«59) 
for  M-(co')  reduces  to  the  free  case^  expression  ['J .hi) ,    when 
the  potential  goes  to  zero. 

Referring  back  to  the  variational  principle  (7.20) 
we  see  that  u(p-,)  =  0  when  the  potential  is  identically  zero. 
When  u(p.)  =  0,  the  f(|p-,l)  — >  1.   The  terms  I^  1    ,    I   and 
Ih  must  all  be  zero  and  we  are  left  with  only  the  term  with  I;-, 

We  have  evaluated  explicitly   I.  ,    1^,    I-.  and  lu    for 
the  specific  potential  U(|p-,|)  =  b  e      .   Suppose  b  =  0. 


-a  IP 
We  must  show  that  I-,  ^  ~L^,    1^  and  I,,  are  in  fact  equal  to 


"1^  ^2'    "^3  '^'■"^  "^4 


2 
zero.  Since  I,  is  proportianal  to  b,  and  II  to  b  ,  these 

two  terms  are  obviously  zero  when  b  is  zero.   Consider  the 

trial  function  f(lp-,|).   As  was  stated  previously  the 

-u(  p-i  1  ) 
correct  f  would  be  e      '  ^  which  should  now  be  1. 
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We  had  assumed  that  f(|p-,|)  =  e   ^      .We  see  then  that 
when  ti.-  i^-lential  goes  to  zex^o,  A  must  go  to  zero  or  Y  must 
go  to  infinity.   With  either  of  these  conditions,  I,,  1^, 
I-,  and  I^  all  equal  zero  identically. 

To  c>-u.i.-^ote  the  argument  wa  must  show  that  the  I 
term,  which      :  rs  in  expression  (7o8)j  reduces  to  the 
free  case. 

In  the  general  case  we  chose  trial  eigenfunctions  for 

t-  ■  tion,  i.e. 

-ar 
X(r  )   =   e 

-pr' 
<f(r')   =  e 

it  we 

ik'r 
X(r  )   =  e 

-ik-r' 
4>(r')  =  e 

When  the  potential  goes  to  zero  X  and  <t>  must  approach 
gas  functions.   B-^ause  of  the  spherical  symmetry 
of  the  trial  functions,  they  can  equal  the  free  functions 
only  when  k  equals  0.   In  this  case  both  a  and  p  equal  0. 
In  addition,  they  must  go  to  zero  such  that  a  =  -p. 
If  we  take  the  limit  of  expression  (7.38)  as  a  — >  -p,  and 
then  as  a  — >  0,  we  can  show  that  (7.58)  becomes  equal  to 
the  free  case  for  k  equal  to  zero.   If  we  want  to  consider 
k  very  small  but  not  actually  zero  we  really  should  not 
compare  (7.38)  with  (7.4l).   Instead  of  (7.4l)  we  should 
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look  at  the  s  state  free  gas  case^  i.e.  the  eigenfunctions 

X  and  4>  should  be  the  s  state  components  of  plane  waves 

and  not  just  plane  waves.   However,  if  we  let  a  — >  ik 

and  3  — >  -ik  as  a  and  p  become  small  we  find  that 

expression  (7.58)  reduces  to  (7.4l)  without  having  to  set 

k  equal  to  zero. 

If  we  take  the  limit  of  (7.58)  as  a  — >  -p  and  then 

as  a  — >  0  noting  that   lim  *_(x)  — >  O^,  we  find  that  only 

x->0 
the  last  two  terms  of  (7.58)  contribute  as  all  the  other 

terms  cancel.   These  last  two  terms  arel 

il   h2)     Ufc-  ax  0.  )         q^8i2mp/2_t\^_^2^  .  /^fK+^^s) 

i7.^2j   KCD  ,UX^    OJ  j 3/2   6 r^ , ,5/2 

icD'b.a3_      -icu'tcD  ico'tcD^ 

.  (2Tri  e       -2Fle       )9(-ta)  )-27rie       0(-tia3  )) 


2' 


[0   -a  )^  ^ 


?  2 
-ill  a  x/2m 


-it^p^x/2m 


(a  -p  )^ 
Since  it  is  M-(.X)'  )  we  are  comparing  and  not  il(<^'  ) 

we  must  multiply  the  whole  expression  by  the  normalization 
factor  (a+p)^/8ir.   Also: 


(P^  -  a^)   =   (6  +a)  (p-  a) 


So 


-86- 


U(|  Pil  )->0  1-^^     h-^  v^ 

X  e  + ^ ^  ~    [l-i(8W^   )  ]e  ]  . 

(3-Kx)^(B-a)^   2P  /    ^m 

Cs         ''■         ''        '    -^p)      terms   we   have 

(7.^4)    ,  11.     .(.. )  =  afsiM.'^'lAl.  „p  J  ,.   ,.  ^^,375^ 

u|p-^|->0  i^/^  h"   yif  d         X        ^    (cD^+cOg)^ 

•  [27ri  e  -2Tri   e  e(-^cD   )-2Trie  9(-*^^p)] 


'-i^'"5<°^^'^^ 


-2;  -—       -\  -2' 


-iii^a^x/2m 


-iii^p^x/2m 


Now  let  p  — >  -a. 


lim  ap  (- -^—     +   -^—  ) 


2 


Also,  as  a^p  — >  0,  since  *(x)  _j,q  =  —m  •   x 


the  error  functions  cancel. 

Letting  a  — >  ik, 

,  -iiL^a^x/2m  -ih^p^x/2m       itiVx/2m 

and        e  =   e  =   e 


4a^   4k2 


Therefore 
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(7.45)  lim  tx(a3>)=ii S/      5     ^ 

potential-^  i"^^      h     /ir 

-i(ai^+a)   )E  2     , 

dto^   dcD^   e  ^      ^     ^      iii^k  x/2m 

.    [2Tri   e"-""  ^"^2   _   2/ri   e  ^   0(-ixa)  )-27ri   e  ^9(-hco  )] 


2' 


where   x  =  ax^cu^/CcD-j^+a)^) . 

We  see  that  this  equals  ^{^'  )f;^QQ' 


8a .   The  Remaining  Integrations  for  the  Uniform  Gas  Are 
Approximated  by  the  Method  of  Steepest  Descent . 

In  the  general  case,  the  ox,    and  co  integrations  could 
not  be  done  analytically  in  closed  form  and  an  approximation 
method  had  to  be  developed.   For  this  purpose,  we  studied 
these  integrations  for  the  free  case.   Although  the  free  gas 
case  can  be  calculated  exactly  by  other  methods,  we  wanted 
to  use  the  same  methods  here  as  we  would  use  in  the  more 
general  case.   Since  we  knew  the  answer,  the  validity  of 
any  approximations  could  then  be  checked.   It  should  be 
remarked  here  that  (7.4l)   |J.^^^p(a)')  is  exact.   Therefore 
the  only  approximation  made  in  the  free  gas  case  is  what  will 
be  involved  in  the  co  and  00  integrations. 

To  facilitate  the  calculation,  we  made  the  following 
transformation". 
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O)^ 


co^ 


9 


Oi^X 


^■]_+aj^ 


=  X  ,   ca,  = 


03^  0) 

-|  -  m  Y  ^   Jacobian  = ~  ,  03+0^= 


It  is  important  to  remember  that  the  paths  for  oll  and  co  were 


r 


^ 


rx 


"\ 


00^ 


w  paths  will  be 


X 


(D^ 


"I. 


With  this  transformation  the  expression  (7,4l)  becomes 
(7.46)  n    (CO. )  -  a!8(2jil.^liV!. 


/  dx  J  dco„  - 


(a>g-x)l/2     -iEj,(-!-)      1 


c^ 


^■-      e 


,  2x  2 
k  ti  X 

2m 


[2(ri   e 


iCD'tlOD- 


-ioo'iioo^ 


ioo'txco^ 


-  2/ri  e 


e(-iiCD2)-2fl-i   e  eC-ixcOg)] 


Consider  now  the  integral 

.1/2 


(7.^7)  J  dx  J  dco^ 


('^o-^)' 


"I 


-.(^  -,)  -^^^^ 


2  e   2m   g(_^   j 


tico' 


where        n  =  —p-   .   (n  is  not  necessarily  an  integer.) 
This  corresponds  to  the  middle  term  of  the  factor  in 


-89- 


brackets  in  (7.46).   The  other  terms  are  so  similar  that 
they  can  he  evaluated  by  the  same  methods  as  this  one. 

The  factor  {ui^-xy^     necessitates  a  branch  cut. 
The  paths  for  x  and  co  should  not  cross  any  branch  cuts. 
Therefore,  we  must  remember  during  the  calculation  to 
choose  the  branch  cuts  in  the  following  way! 


X 


cutX- 


f^ path 


^"^  path 


cu 


03, 


where  we  leave  the  angle  5  undefined  for  the  present. 

¥e  will  now  try  to  approximate  the  x  and  ojp 
integrations  by  the  method  of  steepest  descent  through 
the  saddle  point  keeping  only  leading  terms. 

Making  use  of  standard  techniques  the  method  of 
steepest  descent  yields  the  following  formulas  I 

'a 


(7.48)       [  e=^^^(^)    dt 


zh(t^)      /TTa 


'0 


0 


OD 


J 


2n\  ^('"+  1^     .l>,n 


'2' 


where   a^ 


n 


J^  Z$  '  iT^  1""'  '    ^^^^  -Mt,)-"Ht)  , 


t^  is  the  saddle  point. 


Using 


h(t„)-h(t) 


(t-t^) 


0' 


P 


P 


where  the  A   are  obtained  from  a  Taylor  expansion  of  h(t). 
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we  have  -, 


a 


Wlien  evaluated   eKplicitly,    the   first   few      a^     are: 

!2  ..  15  ;^2  -3  .  2  ;^  ^-1/2 
a^  "  ^  ^1^0    2  ^2-^0 

^4   3.7.9.11  h     -6  tll^l^'i^l^ 


a 


0     2'  '5 


7  ,   "1  "0      ~  2"^ 


+  ^A-''-  +2A^A3)  -|a-\  . 

This  method  is  now  applied  to  equation  (7.^7)* 
We  must  bear  in  mind  that  the  path  of  integration  is  now 
the  one  of  steepest  descent  and  not  the  original  path. 
For  this  to  be  valid  the  deformation  of  the  contour  should 
not  change  the  value  of  the  integral.   If  this  is  not  the 
case,  than  it  is  necessary  to  add  corrections  to  the  path 
of  steepest  descent.  This  point  will  be  discussed  more 
clearly  and  in  more  detail  when  we  know  what  the  steepest 
path  is. 

In  equation  (7.47)  we  have  9(-iia)  )  where  9  is  the 
step  function.   Since  the  step  function  is  defined  only  on 
the  real  axis,  before  proceeding,  we  must  by  analytic 
continuation  define  a  function  which  is  defined  throughout 
the  complex  plane  and  which  reduces  to  the  step  function  on 
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tha  real  axis.   Such  a  function  is! 

9(z)  =  -  arg  (-z),  so 
(7.50) 

©(-iiCB^)  =  J   arg  (ticu^)   =  ^  Imag  (Log  (ho^^))  • 

For  this  function  to  be  uniquely  defined^,  we  must  make 
a  branch  cut.   This  is  done  in  the  manner  already  described. 
Since  this  function  is  not  analytic,  the  choice  of  the 
continued  function  is  not  unique.   However,  we  have  not  made 
use  of  any  specific  representation  for  the  step  function  so 
far  in  the  problem,  i.e.,  we  have  used  only  the  properties 
on  the  real  axis  and  so  any  representation  will  serve  us  now. 

We  must  now  find  the  saddle  point  for  equation  (7.47). 
This  involves  finding  the  saddle  point  in  a  four  dimensional 
space  whereas  standard  methods  usually  treat  only  a  two 
dimensional  space.   Therefore,  this  point  will  be  discussed 
in  detail. 

One  method  of  finding  the  saddle  point  in  four  dimensions 
is  to  treat  two  dimensions  at  a  time  in  a  self-consistent  way. 
This  was  tried  and  led  to  very  bad  agreement  with  the  exact 
answer.   It  is  obviously  more  correct  to  take  the  four 
dimensions  at  once,  but  it  is  also  more  difficult.   However, 
we  worked  it  out  in  the  following  way:   We  rewrite  equation 
(7.^7)  as  2  P  p 


(7.51)  J  dx  J   dcog  e 


©(-iia)^). 
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We  know  that  the  plasma  frequanay  is  given  in  the 
limit  as  k  -^  0.   Therefore  it  is  correct  to  assume  that 

^  ^    «  1  . 


Ep  2m 


When  we  apply  these  sa-     ■  .  :     the  general 
case,  we  will  use  a  similar  condition: 

2.2 
a  n  Xq 

(7.52a)      ^-^^ «  1 

"^0 


^^' V^  "'"^o' 


and 


|(a+p)^ii^x  I 
(7.52b)      2 ""  1  • 

where  x^   and  co^  are  the  saddle  points. 

Aside  from  the  similarity  to  the  free  gas  case  we  have 
no  a  priori  method  of  knowing  that  these  assumptions 
are  correct.   However,  at  the  end,  when  the  numerical  values 
are  known,  these  assumptions  turn  out  to  be  very,  very  good. 
They  are  therefore  justified  by  self-consistency. 

The  saddle  point  is  found  by  taking  the  derivatives  of 
the  exponent  with  respect  to  x  and  to  fi)  ,  setting  these  equal 
to  zero,  and  solving  the  resulting  two  simultaneous  equations, 
For  maximum  accuracy,  we  should  include  the  whole  integrand 
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in  the  exponent,  for  example,  writing  9{-h.(X)  )    as 

In  9(-tcD  ) 
e  .   The  inclusion  of  all  the  factors  in  the 

exponent  was  too  cumbersome  to  handle.   However,  with  the 

exception  of  the  9{~h.cii   ),    all  the  other  terms  were  included. 

The  error  introduced  by  not  putting  0(-tiCO  )  into  the 

exponent  is  discussed  together  with  the  error  introduced 

by  keeping  only  the  leading  terms  in  Section  8b  of  this 

chapter. 

Referring  to   the   exponent   as    E(cd  ,x)    we   have   from 

equation    (7.5I) 

2 
(7.53)      E(cD2,x)    =    -iEp(co2+x   +  —^  +nai^)    +  ^  ln{(^,^-x)-3   In  oj^ 

2 
P    ?  CO  2 

k   11  2  X 

for  -Fr-n-     <<   1,    and  where  we  have  written  — — -  as   co^+x+  th— — r 

n  ^^     ,    ^E    .  i     F       rp    +       2x  ^        _       J_ 

0  "-=   375-  +  3-x"    •       -^    ^F    ^2    +  5^2"^"  +''    ''       -0)2 

„    Se   .     .  „   ,,  ,   2x   ,     X     \        1 
Letting  A  =  — we  have 

QJp  —X 

-  i  Ep  (2  +  2A  +  n)   =  ^ 

-  i  Ej.  (1  +  A)   =   ^ 
These  two  equations  have  the  solution 

(^•5'*)     ^ "  J  - 1 .  -2„  =    rii^  '  ><o "  r;-A  =  nh^ 

F 
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Using  this  l:.  .. ._.  ^..-t^  the  method  of  steepest  descent 
is       eparately  f  :■        ■:  in  order  to  ca" 
coefficients  in  the  expansion.   Also,  in  the  x     ,  the 
above  Xq  will  be  improved  by  including  the  contribution  from 
e      ^        which  was  neglected  in  ths  above  derivation  of  x„. 

From  equation  (7.5I),  the  terms  in  the  exponent  involving 
only  CO  are: 

h{u^^)    -  -i  Ep(a)2  +  x  +  j^-^  +   na)^)  +  ^  In  (^2"^)  -  3  In  co^  . 

^2 
We  evaluate  — ^  h(co  )   substituting  the  above  values  for 


^o  •.   h"(a>      ) .      Using 


VAq 


> 


a         A6 


V  E„n  /n2-l8n+128 

2  F     ' 


the   leading  term  for 


„     h(co  )  h(a)     ) 

(7.55)  e        "^     doi     =  e  V^  an        • 

As  was  already  mentioned  we  want  to  improve  x^-.  to 

2 
include  the  effect  of  the  k  term.   Therefore  we  do  not 

use  x^  and  co_  when  evaluating  h(co  ).   Instead,  we  do 

the  following.   We  had 

.     X  ^2^ 


so 


oi~-x   '  ^  ~  1+K 
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CL). 


^  .n   -  ^iil±Zl 


"^2-^  =  n:+x-    '  ^2  =  — T 


A  was  calculated  to  be  ■^■  ~1,    and  using  the  value  of  A^ 
we  substitute  into  h(cDp) 

^2  =  ""^    '  "^2-^  =  YT^   • 

We  find  for  the  leading  term: 

(7.56)       /rao       ^  ^'^''^^  E^' 


2   /tt  -1   x-^n-^ 


where  N  =  ^  — ^   . 

As  before^  taking  all  the  x  terms  into  the  exponent 
except  for  the  factor  0(-xii  —  jtj-)  ;.  the  exponent  is 

(7.57)         h(x)  =  -  igx  -  ^  log  X  ,    g  -  EpN  -  -2-£-  , 
h»(x)  =  -  ig  -  ^  ,   ^^'(xq)  =0,   ^0  "  2I  • 


If  we  let  k  — >  0  then 

^"F  ^  n^     -^ 
which  is  the  same  as  the  first  Xq  we  found.   Similarly 


(7-59) 

^  IT  ■?  On 

as  k  — >  0, 
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^2q      -      ^O^nV           ^^^ 

2i 

^F  5-  n-4 

V 

which  :      .;ame  as  the  first  a>^  ,  i.e.  expressions  (7.5^). 

Therefore,  v,   '  . ..  —  ^    •  ^ 

include  the  effect  a. 

The  leo  .  .'.^    i  k  integration  is 

N  r  A/  ^   h(x)  ,     ^1       N   ^(^0^   r^ 

.   )         J  <J)(x)  e  ^  ^  dx  =  <f(xQ)  e      /ir  a^  , 

where         a^  =  -"'---  ,   ^o  '"  "  ^   ^^^0^  ®^'^  ^0  "  g"  ' 
The  111  _:  .^    ;       ^       (T-'^T)  becomes: 

V  5/2 


2 e    9{-   ^n  ^) 


\_   5/2         ^  n-^  2s' 


-to 


Referring  back  to  the  expression  (7.46)  for  M.(a)' ) 

free 

we  note  that  the  last  term  is  the  same  as  this  term  with 
n  — >  -n.  The  first  term  is  also  the  same  if  n  — >  -n  and 
i--    ---^.J    ^  1. 

Therefore   the   leading  term  for  l-'-('^' )f  j^gg    is: 

il  62)    11(0)'  )  =  a_°l2m) ^j       e_         2 

U.D^;    M-^u)   ;^^g^  ^5/2  ^6    ^  j^2  25 

.[2Tri   K^   g|/2   -2.iK^   gj/^   0(-  |n_  |J^)    _^^^  ^^  ^3/2 

Q/      hn        5  i\ 1 

,  2j^2  j^2^2 

where  g3_  =   E^N^   -  -^^  ,        g^  =   E^Ng   -  -g^T  ' 


Ejpn^  y^-l6n+128  E^n^  v^^+l6n+128 
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8ti .   The  Error  Introduced  by  this  Approximation. 

We  must  now  discuss  the  error  Involved  In  keeping  only 
the  leading  term  and  also  from  not  Including  the  step 
function  In  the  exponent  when  evaluating  the  saddle  point. 

Using  standard  methods^  the  first  few  terms  in  the 


asymptotic  approximation  of   dco  ^(cDp)  e 


2' 


are . 


(7.63)   <f(^o  )  e 


20' 


0 


a   /^  n  ^  (2-1)1  J   2  ^  3:1 

a^  -A  [1  +  ^-^-   ^  +  V  ~  +• 


'4 


0 


'0 


+  *'(a)2Q)  e     0  ^   t  _^   [2-11^ 


+ 


*''(a^  )   h(a)2Q)       a^(2-l)ll 
+ ^-*  e      /r   [ 5 + 


where  *  '  (co   ) 
"^0 


*"(^o  ) 


_d_ 

dCLi^ 


Iraag  (Log  (iio)  )) 


TT 


0^2=^2 


d^   1 


2 


0 


da3^ 


TT 


-  Iraag  Log  (iioj  ), 


^2=^"2 


0 


Since  the  Imag  (Log  (z))  is  not  analytic,  we  must  say 
what  we  mean  by  the  derivative.   We  can  consider  this  function 
as  a  real  function  of  two  variables  and  take  the  directional 
derivative  tangent  to  the  path  of  integration.   When  we  look 
specifically  at  the  directions  involved  we  find  that  the 
result  is  equivalent  to  taking  the  imaginary  pert  of  the 
derivative  of  the  analytic  function  Log  (z).   In  other  words, 
even  though  it  is  usually  not  valid  to  interchange  Im  and  -^— , 
in  our  case  we  found  that  the  two  operations  did  commute. 
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It  is  quite  straightforward  to  e'/aluate  the  coefficients 
a.  and  the  derivativea.   We  find  that  for  s.-nall  n,      i.e. 
n  less  tl       .    !  ratio  of  the  terms  witxh  *'  and  <J)''  to 
1         term  with  <t>  is  0.12  independei.   ^   f  n. 
H  ■"  '-^r,  for  larger  values  of  n,  we  find  that  this  ratio 
becomes  dependent  on  n.   In  fact,  it  decreases  as  =—  . 
We  defined  :.     .e  —^-   .   F^         _:as  cu' ,  the  plasma 
f "  — " -ncy,  is  proportional  to  the  square  root  of  the  density. 
The!-     J  n  increa      ■.  the  sy:      ecomes  n, 
the  above  ratio  decreases.   For  a  f^  jas  at 

rr.  ■^'^""'ic  densities  n  is  about  2  and  the  ratio  is  about  0,02. 
For  a  free  electron  gas  at  ■    '  :  densities  the  ratio 
I        .    .   If  we  use  the  value  of  n  for  an  atom  that 
we  calculate  at  t'- -   id,  this  ratio  becomes  'vO.  00002. 
This  ratio  repre.       .  j  error  introduced  by  not  including 
<t»(cu  )  in  the  calculation  of  the  saddle  point.   We  see  that 
for  very  small  densities  the  maxiraxom  error  is  12  percent 
and  for  very  dense  systems  the  error  becomes  negligible. 
Since  we  are  interested  primarily  in  very  dense  systems, 

this  is  a  very  good  approximation. 

ap  a^. 
We  find  that  the  ratios —  ,  —  do  not  depend  on  n, 

^0  ^0 
nor  any  other  parameter  in  the  problem.   These  numerical 

coefficients  decrease  fast  enough  by  themselves  to  warrant 

taking  only  the  first  term.   However  we  find  that 

1^1  =  0.04    and    iS^Iii  |!2|  =  0.02  . 
^0  0 
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Although  this  is  small  enough  to  justify  keeping  only 
the  first  term,  it  limits  the  accuracy  we  can  expect  from 
our  final  answer. 

Repeating  the  same  procedure  for  the  x  integration, 
we  again  find  that  the  ratios  of  the  terms  in  the  asymptotic 
series  to  the  first  term  depend  only  on  the  numerical 
coefficients  which  fall  off  quite  rapidly. 

I®2i  1  |^4|  1 


IbqI  20     >  la     I         1000      ' 

and 

(2-1)11    1^21  1 

"2  IFqI    -  ZTo     =   0-025        . 

The  total  error  is  less  than  or  equal  to  the  sum  of  the 

errors  from  the  cOp  part  and  the  x  part.  Therefore,  the 

error  introduced  by  keeping  only  leading  terms  should  "be 
less  than  4 . 5  /b. 

8 c .   Comparison  of  the  _Pa_th  of  Steepest  _De_SG en t  _w_l t_h  the 
Correct  Path  of  Integration. 

Since  we  are  keeping  only  the  leading  term,  the 
direction  of  the  path  of  steepest  descent  is  just  the 
direction  of  steepest  descent  in  the  region  of  the  saddle, 

The  condition  for  this  path  is 

o 

(a)  h"(co   )  (co  -  CD   )   :=  negative  real  number 

^0   '^    "^0 

(b)  h"(x„)  (x-x^)      ■=  negative  real  number. 

Upon  evaluating  h"(<X)   )  and  h"(xp)  we  find  that  they 
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2 
ar>^   bnr.h  negati-.  -    ■_::-_   ..,.::.bers.      We  v;rite    (od  -co      )      a: 

0 

0*^0  and   h"  '         ''    as      R|e      .      From   equ- •■ 


7r  +  29-Tr,  9-0, 

so  that  the  path  of  steepest  descent  in  the  to  plane  is 
just  a  horizontal  line  through  the  saddle  point.   Wj  find 
the  same  result  for  the  k  plane. 

The  original  paths  of  integration  were  also  horizontal, 


— > 


^2 


In  addition,  they  lie  above  the  real  axis.   This  is 
important  because  there  is  a  singularity  at  the  origin. 
We  have  found  that  the  saddle  points  were 


^^   -v  -FT^   ,  x^  = 


1  8(^  -1) 


Consider  just  the  oo  plane  as  similar  remarks  will 
hold  for  the  x  plane.   If  n  is  real  and  positive,  then  a)_ 
lies  on  the  positive  imaginary  axis.   In  this  case,  the 
path  of  steepest  descent  coincides  perfectly  with  the 
original  path  and  no  correction  is  needed.   However,  in 
our  problem,  we  do  have  terms  with  n  a  negative  real  number, 
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In  such  cases  co,,   lies  on  the  negative  real  axis  and  the 
steepest  path  is  then'. 


We  note  that  this  path  is  on  tha  wrong  side  of  the 
singularity.   Also,  it  is  not  permissible  to  cross  the  branch 
cut.   Therefore  J  we  must  add  correction  terms  to  this  path. 


Choosing  the  angle  5  infinitesimally  small,  the 
contributions  for  J-,  and  J^  are  just  that  from  the  path  of 
steepest  descent.   However,  the  contribution  from  Jp  must 
be  added  on  to  the  value  for  the  integral  previously  derived 

The  evaluation  of  the  contribution  from  Jp  is  just  as 
hard  as  doing  the  original  Integral.  However,  the  problem 
was  avoided  in  the  following  way.  In  the  free  gas  case  we 
know  the  correct  answer  and  we  say! 
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correct  answer  -  expr  :  •  ■  ' -n  (7.62)  +. correction  from  J^. 

We  ev" '■:-'■»  the  expression  (7.62)  in  the  limit  k  — >  0, 

•le  exact  answer  in  the  same  limit. 


2   _  _^_   ^fL-P-l^  _  2_^_^_  1_^"I. 


"p 


Ep 


^ ^-^-  ==  ^2717^^7^^"  '      ^   ^'  ^^^  density; 

-/5 


jxpression  for  (7.62)  for  m.((jo')  as  k  — >  0  has  two 

2   2 
terms.   One  becomes  ] ,„a.--„,^-|_  ^^  ^n  Ai  .   The  exact 

answ-ar  in  this   1 ;  '  .    .   econd  term  is 

singular  in  lA  •   I-      exact  answer,  there  is  also  an 

.  2 
apparent  singularity  in  lA  ,   But  as  k  — >  0  the  numerical 

coefficients  in  the  exact  case  are  such  that  the  terms  with 

2 
lA  all  cancel. 

Therefore,  the  expression  (7.62)  differs  from  the 

exact  answer  as  k  — >  0  not  in  any  physical  constants, 

but  only  in  numerical  coefficients.   Therefore,  the  form 

of  the  correction  due  to  Jp  must  be  the  same  as  the  form 

of  M-('-j' )  in  (7.62),  but  with  numerical  coefficients  such 

that  (7.62)  +  Jo  =  exact  answer.   After  doing  the  arithmetic, 

this  numerical  factor  is 


[-  1  +  c  i]  , 
where 

-105- 


^        12.5   F 

3(-n^-4)2  ,        ^)2 

F  =    [ 1 _—  +  ^     P    --^-=-    ] 

n      yii2+i5n  +128  n      /n2-l6n   +128 

.=  ^-i|/^        ,  p  .   density; 


for  metallic  densities   c  =   O.O2183 
for  atomic  densities    c  =  O.O2836. 


The  contribution  due  to  the  J„  term  is  just 


(7.64)         [-  1  +  c  i]  •  |i((u'  ) 


from  (7.62)  * 


9 .   The  Integration  Completed  for  the  General  Case. 

¥e  now  refer  back  to  the  expressions  (7.59);>  (7.56), 
and  (7.38).   It  is  desired  to  complete  the  calculation  of 
il(a)'  )  by  performing  the  integrals  over  ecu  and  -Dp.   In 
(7.36),  I-,  and  Ip  are  independent  of  cu  and  ^a3    Define 


a  z^  and  z^   so  that 


I3  --=  (^1  +^2 ^^5 


(a^  +  ^i|) 


where  z_  and  Zl  are  independent  of  ecu  and  '^13  . 

We  repeat  the  same  methods  now  as  we  used  for  the  free 

gas  case.   The  calculation  falls  into  two  parts.   The 
first  is  the  method  of  steepest  descents.   Apart  from 
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sligh^''  -  "lore  difficult  - '  "'  ■■  .  the  procedure  is  exactly 
the  sa.-ne  as  before  and  needs  no  further  discussion.   As 
was  men;     i  in  the  free  gas  case  though,  we  did  assume 
here  that 


(«)  E^rL-   -   1 


(7.52) 

,  apart  from  the  e(-iiu)^)  fu:.  ..,  it  was  p^  .e 
'^  -  ■  -elude  th'  ■  ■^ire  function  in  the  saddle  point  calcula- 
ti  .      J  saddle  point  differed  for  the 

different      ,  .)   To  include  the  error  function, 

•-"      ^  '  '^    '     '     approximations  (a)  and  (b)   and 
.0  approximation 

<J>(x)  -w  -^-  e"-"^  [X  +  ^-  +   . . .  ] 

which  is  valid  for  small  x. 

The  second  part  concerns  the  correction  to  the  path 
of  steepest  descent  when  that  path  was  on  the  wrong  side 
of  the  singularity.  We  laiow  what  this  correction  had  to 
be  in  the  free  gas  case.  Now,  when  we  compare  the  integrands 
in  fi  CD' )  in  the  general  case  with  equation  {J.hj)  we  see  that 
the  different  terms  in  ii(<D'  )  differ  from  equation  (7.47) 
only  in  the  following  ways; 
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(1)  k  is  replaced  by  either  -a  or  -(a+p) 

(2)  sometimes  there  are  positive  integer  powers  of  x  . 
as  extra  factors 

(3)  sometimes  there  is  an  extra  factor  of  y^T 

{k)      the  terms  with  z  have  the  additional  factor  (co  +cd  ). 

We  should  remark  here  that  the  terms  with  Zk  do  not 
have  a  singularity  in  oOp  or  x  and  so  we  need  not  consider 
any  correction  to  these  terms. 

We  can  rewrite  equation  (T-'^T)  3S 


0),- 


(7.65)   /  dx  ^ 


(^  _x)l/2   -iEpCn^^-  +n.a)J  ikVx 


2 


"•"2  -■  0^ 


2 


2m 


eyx 
lim  e      (-iitD  )  . 
e->0 


The  correction  was  evaluated  over  a  path  Jp  in  the 
X  and  o)  planes. 


On  the  path  Jg,  all  the  terms  of  the  integrand  of 
(7.65)  are  bounded  as  is  the  path  length  on  :S ^.      Therefore, 
equation  (7.65)  is  uniformly  convergent.   Uniform  convergence 
is  sufficient  to  insure  that  it  is  exactly  correct  to 


-106- 


interc:.  :,_  limits ,  ir.   _■  '.ions,  and ,  j  i.    l 

To  find  the  correction  to  equation  (7.65)  due  to  J^,   we 

re-evaluate  the  free  gas  case  with  the  e.Ktr- 

lim  e   '    in  the  u.   _ 

We  find  that  the  correction  term  for  equation  (7.65) 


now  becomes 


fir  1       Lp^  +f^J 
(7.66)      lim      [-l+ci]7r   K  /^  ^     --§ 

e  ->0  ^  -^         r5i 


Si      "^S[||  +a£]  e/H 
[■gi  +Y£]e    ^S     e  V^S 


#  ^  -1^'^ 


iin   /5i 


20 


g  =  gj_ 


g  =  §2 


if  n  >  0  ;         Sif  n  <  0  . 

Referring  to  equation  (7.65)  we  can  evaluate  the 
corrections  to  each  of  the  terms  in  the  general  case  by 
the  following  procedures. 

(1)  replace  k  by  -a  or  -(a+p)   as  the  case  may  be 

(2)  if  there  are  extra  powers  of  x,    differentiate 

2  2 
equation  (7.66)  by  d/d[k  h.    /2m]  for  each  power  in  x. 

This  is  equivalent  to  d/dg  in  equation  (7.66) . 
(5)   if  there  is  a  factor  yx~  differentiate  (7.66)  with 

respect  to  e  before  takimg  lim  . 
(4)    -iEp(--f--  +no^^)  -iEplcOj^-Kiig)   -io^'co^ 
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So  for  a  factor  {a^i  +(X)   ) ,    differentiate  equation  (7.66) 
with  respect  to  £„. 

This  procedure  is  quite  lengthy^  but  quite  straight- 
forward and  not  at  all  difficult.   However,  it  has  the 
obvious  disadvantage  that  taking  many  derivatives  generates 
many  terms.   Therefore  the  final  result  for  M-(co' )  in  the 
general  case  is  very  cumbersome. 

After  carrying  out  all  the  procedures  indicated  in 
the  last  few  pages  we  find  for  p.(a)'): 

(7.67)  fi(co')  =  [2T(-n)  -  2U(+n)  -  2U(-n)]  , 

where  T(n)  =  U(n)  with  all  the  0  functions  equal  to  1 
and  all  the  derivatives  of  9   equal  to  0. 

2^5  Y^  m^/^  q2p  iap[I^+l2] 


(7.68)        U(n) 


h^  yW 


) 


.Vci  v^Ke5/2g3/2e(^|_ 
2^^ir^{a+^)c  i  K  ^  e^^^  g^^^  d{g^) 
2iiV   G    i    K  y^  e5/2    (EpN)^/2    ©(E^N) 


t\^K  gj/^   2^/2   e5/2  ^(^^  y57r5mg£)    e(g^) 

2^   (3      4m   25 
hV(a+p)K  g^/2  25/2  ^5/2  ;^((ct+p)^  v/5A6m"g2  ^(gg) 

2^  (Tp^  4m  25 
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hV'^^    K  e^/^   EpK    z    ©(E^W)     (3a+f3) 


F 


+ 


^/^(a+p)h5   Ky^EpNe^   e(-^--     ^i^ 


TF^"p^"l6  m^"^^' 


EyN' 


(a+B)2iiV   y^K  5   &\^^   e5/2(_i-ci)9(g2) 
^  "a'"2l°  p5   8m2 


+ 


ic(a+p)2iiV   ^K  3e5/2    (e^n)^/^   9(EpN) 
~  25  a  2^°  p^   8m2 

(a+p)  Vh^   K   y^  g^/2   e5/2    9^^(g    )     (-1-ci) 

25  a  2-'-°  p^  4m2 
ci(a+p)2tV  K    v^   (EpN)5/2   ^5/2    0g(EpN) 

25  a   2"^°  p^  4m^ 

(a+3)^  K  7r^/2  ^5  e5/2   ^   ©(e^n) 
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(a+p)^  K  Tr5/2  ^,5  ^^^  ^5/^   ^   ^g^V^ 


"a"2l0  p5  J,  ^57^25 


(a+p)^    K   ^5/2^5   yrf  e5/2[l(E^N)^/2    9    (en)    +(£,4^)^/2.   ^^^^ 


"a  710  p5"4  m572   25 


?g_Z_-i  . 


15^5/2q2^p^3i  V^  z 
+ 


h^  )/ir 
-iii^-fr^L     v^  e 


3/2   ^1/2   g.iin       3i^ 


=1 


2        p     m  3 
.2.  '2        7:7  ^5/2  K  g^/2        g3/2   ^  Q 


16   p      m   25 


E    2S 


e: 


h^yS,    ^K   e5/2   g3/2    0^    (       ) 


1 


2^   B      m  25 
1  L  k=^  ir^  v^  gl/2  e3/2  /- 


r 


^n-T  2g.  ^ 


2^2, 


a   2^   p^  m  25 


~io     A 
ct   2^^  3      m  3 

5/2  K  g^/2      g5/2n 


Ep   2- 


E^  2 


F 


2^2^,/o«5/2  K(E^N)5/2         (EM)5/2n 


ii^TT^cv^  e 


Q 


a  2     p-^  m  25 


^F    2 


5 


^ 


2|      (EpN)l/2p]0(E^^) 


^V(a+P)c    >^K    e5/2[-gV2    g^    (       )    _    (e  n)5/2   q      (E„W)] 


+ 


_      ___ 
2^   a  p      m   5 


Ep^-F^ 
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f,      '-\ 


T        /r^  ^      /o 


L   \. 


aV^'7''m7 


+ 


ihVo 


ihV/2 


^/2 


•O 


-  i(  (-«)>> /irk- )«(^-K- 2-1^) 
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CliAPTER  VIII 
Description  of  Numerical  Methods 

Tlie  numerical  work  was  done  only  for  method  I.   In 
this  case,  we  have  a  rather  complicated  expression  for 
M-((jci'  )  which  must  he  stationary  with  respect  to  the 
parameters  a,  ^,    and  also  with  respect  to  those  parameters 
appearing  in  I,,  Ip,  I  and  I l . 

In  order  to  do  the  numerical  work  it  was  necessary  to 
choose  a  specific  form  for  the  external  potential.   We 
decided  that  an  atom  would  be  an  interesting  example  of  an 
inhotnogeneous  electron  gas  to  study.  The  external  potential 
for  an  atom  would  be  the  Coulomb  potential.   However,  when 
the  problem  was  formulated,  we  did  not  actually  calculate 
the  self-consistent  correction  to  the  potential  as  we  should 
have  dona  because  of  the  difficulty  involved.   However,  we 
know  that  if  we  had  made  this  correction,  the  external 
potential  would  be  the  Hartree  potential,  which  would  be 
a  screened  Coulomb  potential.   Therefore,  it  was  decided 
to  choose 

V(r)   =  b  e"^^  . 

For  typical  atomic  parameters  we  chose  those  for 
Krypton.   For  Krypton  the  Fermi  level  (the  energy  of  the 

highest  bound  state)  is  -  22.4  x  lo""^^  ergs.   We  let 

?   -R 
b  =  Zq  AO   .   For  Krypton  Z  =  36  so  b  is  II5.27  EL,. 

o 
We  first  chose  a  to  be  10  and  then  for  a  confidence  factor 

Q 

we  let  a  be  10  . 
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_  Q  V> 

For  a  potential  of  the  form  be    ,    we  have  already 
computed  I^,  I^,  I^  and  I^  (see  (7.52),  (7-53),  (7-5^) 
and  (7.55)).   These  terms  Involve  the  parajneters  A  and  7. 
These  expressions  I, ,    Ip,  I^  and  lu   must  be  substituted 
into  lJ.(co'  )  . 

|-l(co'  )  must  now  be  stationary  with  respect  to  a,  ^,    A 
and  y.      Therefore,  we  took  the  partial  derivatives  of  M-(oo' ) 
with  respect  to  each  of  these  parameters  and  required  these 
four  derivatives  to  be  zero.   Since  the  derivatives  with 
respect  to  a  and  p  were  rather  lengthy  and  so  error  prone, 
they  were  checked  numerically.  The  derivatives  of  I-l(co' ) 
with  respect  to  a  and  p  at  a  point  (a^,B^)  were  evaluated 
numerically  by  a  difference  method.   These  numbers  were 
then  checked  with  the  analytical  expressions  for  these 
derivatives  evaluated  at  the  same  point  (qq^Pq).   In  this 
way  we  could  be  sure  that  the  analytical  expressions  for 
these  derivatives  were  correct  and  also  that  they  were 
programmed  correctly.   This  same  method  also  checked  that 
|J.(co' )  was  correctly  programmed. 

The  condition  for  a  resonance  was  Real  (ti.(a)'))  =  1. 
The  imaginary  part  of  |J-(a)'  )  was  a  measure  of  the  damping. 
Since  we  were  interested  in  only  observable  resonances  we 
required  zero  damping  --  i.e.  that  Imag  ((i.(6o'))  =  0  for  a 
real  cd'  . 

We  therefore  had  to  solve  simultaneously  the  following 
six  equations: 


■iia- 


RFl  =  Real  M.(a)'  )  =  1 
IFl  =  Imag  M.(co'  )  =  0 


(8.1) 


P2           : 

aa 

=    0 

F3      = 

=    0 

F4             : 

=    0 

F5     -■ 

-    0 

for  the  unknowns  n,  a,  p.  A,  Y  where  n  ~ 


^F 
We  then  defined  a  function  F  to  be 

F  =  (RFl-  1)(RF1-  1)*  +  (lFj_)(lF^f  +  c^F.^?^ 


where  the  c's  are  positive  numbers  and  the 
asterisk  signifies  complex  conjugation. 

We  see  that  a  zero  in  F  corresponds  to  a  solution 
of  all  the  equations.   Furthermore  we  see  from  the  definition 
that  F  _>  0.   Therefore  to  solve  the  system  we  minimized 
F  hoping  to  find  that  at  the  minimum  F  would  be  zero. 
Since,  in  method  II  we  found  that  the  approximations  made 
had  an  error  of  4.5  percent,  we  decided  to  count  as  a  zero 
any  values  of  F  that  were  zero  to  J>   or  more  significant 
figures.   The  c's  were  introduced  as  scaling  factors  so 
that  each  of  the  separate  equations  would  have  equal  weight. 
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If  we  look  at  ths  specific  expressions  for  F4(^) 
and  F5(-v^)   we  see  that  f4  has  a  Y  in  the  denominator 
and  F5  has  ay.   We  found  numerically  that  these  denominators 
dominated  these  equations  and  we  were  getting  solutions  for 
7  — >  00  ,   We  have  shown  before  (see  page  85)  that  in  this 
case  M-((X)' )  reduces  to  the  free  gas  value.   To  avoid  this 
solution  we  therefore  replaced  F-^  and  F5  by  the  numerators 
of  f4  and  F5,  i.e.  we  required 

numerator  f4  =  0  j         numerator  F5  =  0 

This  resulted  in  a  modification  of  F. 

F  =  (RFl-  1)(RF1-  l)""  +  IFilF^*  +  C2F2F2* 

+  c^F^F^*  +  Y^c^F,^F^*  +  Y^c^F^F^*  . 

To  minimize  F  we  used  a  routine  called  BOTM  from 
the  GIMS  Library.   This  routine  had  the  advantage  of  not 
requiring  derivatives  and  also  it  converged  well  even  for 
very  poor  initial  guesses. 

It  must  be  mentioned  that  the  following  restrictions 
on  the  parameters  were  ma  del 

(1)   a^p  occur  in  the  expressions  for  eigenf unctions 
e"^'^  and  e"^"^.   Since  it  is  unphysical  for  these 
functions  to  become  infinite  as  r  goes  to  infinity^ 
we  required  that  the  real  parts  of  a   and  ,3  not  be 
negative. 
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(2)  If  the  real  part  of  A  became  more,  negative  than  -100 
then  the  value  of  F  hecame  too  "•  ■   ■•   for  the  machine 
to  handle.   Since  we  were  looking  for  a  minimum^  we 
knew  that  the  solution  for  A  could  not  lie  in  that 
region.   Therefore  the  real  part  of  A  was  not  allowed 
to  be  less  than  -100. 

(3)  We  were  looking  for  a  positive  frequency. 

Since  E„  was  negative  and  n  -  -pr-  ,  we  required  n  to 
be  negative. 
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CHAPTER  IX 
Results 

Two  distinct  solutions  of  equations  (8.1)  were  found. 
For  the  first  solution  we  determined  the  values  I 

n  =  -  16.3 

a  =  (1.58  X  10^,  -4.87  X  10^)    ^^ 

p  =  (1.18  X  10^,  9.26  X  10°) 

Y  =  (6.47  X  10^,  -1.72  X  10^) 

A  -  (127.5.  -35.4) 

hoD'  =  En   -    22.4  X  10    X  16.3  = 

=  3.67  X  10""^°  ergs 
^  =  54  A 

This  solution  was  very  unstable  with  respect  to  changes 
in  a  or  p  of  0.1  percent. 

This  highly  singular  solution  reminds  us  of  the 
breathing  mode  found  in  a  system  of  harmonic  oscillatoj^s 
which  are  coupled  to  a  common  center  and  to  each  other. 
The  second  solution  was  not  so  sharp.   We  found  the 
following  ranges  in  the  parameters  corresponded  to  a 
solution; 


1'^ 
^   The  two  numbers  in  parentheses  refer  to  the  real  and 

imaginary  parts  of  the  parameters. 
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n  -  -  1.78  X  10^   to  -5.70  X  10^ 

a  (5.28  X  10^  to  1.26  X  10°,  -^1.68  X  10  to  -4.56  x  10^) 

p  (5.27  X  10°  to  5.80  X  10  ,  7.64  x  10   to  8.59  X  10^) 

A  (8.52  X  10^  to  4.0^  X  10^,  -1.21  X  10^  to  -1.27  X  10^) 

y  (1.18  X  10^  to  1.13  X  10°,  1.30  X  10^  to  1.98  X  10^) 

There  is  a  reason  for  the  lack  of  sharpness.   We  found 
that  for  the  above  ranges  in  the  parameters  a,  p,  equations 
Fp  =  F^   i.e.  -c—  =  -^  .   Therefore,,  instead  of  uniqueness 
in  this  region  there  is  a  family  of  solutions  since  the 
problem  is  underdetermined . 

8  R 

tcD'  =  EpH   ranges  from  3.99  x  10"  to  8.29  x  10~   ergs 

and 

Q-n-,^  o  o 

— —        ranges  from  0.24  A     to   0.46  A  . 

These  numbers  would  indicate  a  K  shell  vibration  rather 
than  the  whole  atom  and  the  lack  of  sharpness  would  indicate 
that  the  vibration  of  this  shell  is  perturbed  by  the  other 
shells.   Our  original  ansatz  that  the  solution  to  our 
integral  equation  be  an  s  state  may  be  forcing  an  s  state 
solution  which  is  not  really  correct.   This  shows  up  as  a 
solution  which  is  not  sharp. 
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CHAPTER  X 
1 .  General  Remarks 

We  have  shown  that  two  separate  collective  states  may 

,  o  ,0 

exist  in  an  atom;,  one  at  54  A  and  the  other  between   .24  A 

and   .4o  A.   There  may  be  other  solutions  in  addition. 

Our  solutions  were  found  from  the  zeros  of  l-l^(cu'  )  where 

|J-(cD' )  is  an  eigenvalue  of  equation  (2.53)  corresponding 

to  the  3  state  eigenfunctions  <t)(x')  =  e"^'^  '  and 

^(x)  =  e"'^''^'.   There  are  other  eigenvalues  corresponding 

to  higher  order  eigenfunctions.   It  may  be  that  some  of 

these  other  eigenvalues  will  lead  to  additional  collectxve 

states.  Furthermore^  our  final  solutions  were  obtained  by 

numerical  methods  which  unfortunately  cannot  guarantee 

uniqueness.   Therefore,  a  more  extensive  numerical  search 

might  discover  additional  collective  states  even  for  the 

|J.(cD')  that  we  have  worked  with. 

An  additional  point  needs  to  be  discussed  concerning 

the  validity  of  our  results  with  respect  to  a  real  atom. 

In  our  derivation,  we  considered  a  nonuniform  dense  electron 

gas  perturbed  by  the  interelectron  Coulomb  force.   The 

nonunif ormity  was  due  to  a  self-consistent  external 

potential  which  we  have  shown  to  be  equivalent  to  a  Hartree 

potential.   However,  in  our  application  to  an  atom,  we  did 

not  use  the  exact  Hartree  potential  for  an  atom.   Instead, 

we  approximated  the  density  variation  typical  of  an  atom 
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by  choosing  the  form  of  the  external  potential  to  be 
be    ,    with  the  amplitude  b  and  the  range  a  chosen  to 
approximate  an  atom.   We  have  answered  the  question  as 
to  whether  very  nonuniform  systems  can  still  have  plasma 
oscillations.   In  fact,  we  have  shown  that  an  atomic  system 
should  show  these  effects.   However,  because  an  exact 
atomic  Hartree  potential  was  not  used,  our  results  are 
probably  valid  only  as  an  order  of  magnitude  estimate. 

2 .   Comparison  with  _0t he r  Theoretical  Work  in  the  Subject . 

Unfortunately,  there  is  not  much  other  work  to 
compare  with.   In  fact  to  our  knowledge,  there  does  not 
exist  any  work  besides  ours  which  includes  the  damping 
effects.   However,  we  will  discuss  those  few  papers  that 
there  are  which  deal  with  the  subject  of  collective  states 
in  atoms. 

Several  papers  in  this  subject  have  been  published 

by  W.  Brandt  and  S.  Lundqvist.   In  the  first  of  these 

l4  ^ 

papers,    the  authors  consider  a  response  function  R(q,CD) 

defined  as 




W.  Brandt  and  S.  Lundqvist,  Phys .  Rev.  1^,  21^5  (I963). 
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(10.1)  R(q,a))  =  J  d^x  d^x'  e""-^*  (''"''' ^  ^{x,x^  ,o^)    , 
with  Q  given  by  the  integral  equation 

(10.2)  Q  -  Qq  +  QqVQ 

Q„  is  defined  as  auto-correlation  function 

in  the  Hartree-Fock  picture  and  the  expression  given 

for  Q(^  is  exactly  the  same  as  our  expression  for  V  in 

lowest  order.   (See  equation  (S.-^).)   The  resonances 

of  the  system  correspond  to  the  poles  of  R(q_,a3).   We  note 

at  once  the  correspondence  between  (10. l)  and  (10.2)  and 

our  equations  (2.27)  to  (2.52). 

The  authors  try  to  take  into  account  the  effects  of 
a  density  gradient  when  evaluating  Q^  by  using  a  formula 
due  to  R.  Glauber.   In  so  doing;,  they  find  that  Qq 
depends  upon 

(10.5)  <\[f   -  Pf(x)> 


and   that 


oo 


,i(zt-t^) 


(10.4)  <0  (z)>=  ±  pk  i^'^     J  —--= — -     ^^  • 

±  ^"^    e-^0  t    +   is 

-oo 

This  is  equivalent  to  the  approximation  we  have  for 
the  density  operator.   If  they  had  proceeded  along  these 
lines ;,  they  would  have  done  what  we  have  done.   However^ 
in  solving  the  integral  equation  (10.2)^,  Brandt  and 
Lundqvist  use 
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whi :     the  same  as  if  we  were  to  use  plane  wav-'S  for  our 
fun:  :        i  V'-   They  then  find  that  in  the  long  wave 
length  limit  expr  •''  ''10.4)  gives  the  same  result  as  a  gas 
with  no  density  gradient.   This  is  not  surprising  to  us 
as  we  found  that  our  result  reduces  to  the  uniform  case 
when  we  let  <t>   and  V  approach  plane  waves  and  take  the  ^  -  - 
wavelength  limit. 

I:     /aar  paper^      -ndt  and  Lundqvist  make  u     :' 
the  ■   ^'   '  t  in  the  '  '    quantum  number  limit  the  proper- 
ty a  can  be  st      oy   considering  an  assembly  of 

c-.-i- -  ,.^-- _-  ..u___: _^  _:.  .--_jh  the  coupling  constants 

ean  radii  of  electron  orbitals.   Using 
this  model  they  derive  a      r  relationship  between  the 

total  a_j-._j  moment  ,_  atom  and  an  externally  applied 

field.   From  this  they  determine  the  atomic  polarizability, 
a(a))  in  terms  of  the  single  particle  oscillator  strengths 
f  .  excitation  frequencies  cd  and  coupling  constants  A  . 
The  resonances  of  the  system  occur  at  the  poles  of  a(a)). 
With  all  the  ?v  taken  at  some  mean  X  they  find  the 
disperson  relation 

(10.6)  1  +  A^a(cD)   =   0  . 


-^  W.  Brandt  and  3.  Lundqvistj  J.  Quant.  Spect.  Rad.  Transfer, 
Vol.  4,  pp.  679-689  (1964). 
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This  is  equivalent  to  our  condition  that  the  real  part  of 
the  dielectric  constant  be  zero.   The  solutions  of  (10.6) 
can  correspond  to  single  particle  as  well  as  collective 
effects.   We  found  that  we  needed  a  second  condition  i.e. 
that  the  imaginary  part  of  the  dielectric  constant  be  zero, 
in  order  to  determine  the  collective  states.   By  looking 
at  contributions  to  a(co)  from  the  different  atomic  shells 
and  by  qualitative  arguments,  Brandt  and  Lundqvist  deduce 
that  collective  effects  should  be  in  the  frequency  range 
1  ^    (r— )  ^   Z   where  iiR  =  15. 6  eV.   This  would  put  co  in 


•Ry'   -        -y 


O  ,   o 

the  range  of  10  -  100  A.   Our  first  solution  of  54  A 
agrees  very  well  with  this  result.   Our  second  solution 
does  not.   However,  at  wavelengths  shorter  than  K-shell 
absorption,  the  dipole  approximation  of  Brandt  and 
Lvmdqvist  does  not  suffice.   Therefore,  it  is  understandable 
that  their  method  does  not  yield  a  solution  at  very  short 
wavelengths.   We  should  remark  that  one  of  the  conditions 
Brandt  and  Lundqvist  use  to  determine  which  of  their 
resonances  is  collective  is  that  there  be  "a  sufficiently 
high  gap  on  the  high-frequency  side  of  the  single  particle 

spect3fTJim  of  accommodate  this  new  resonance".   Our  solution 

o         o 
of   .24  A  to   .46  A   satisfies  this  condition  as  it  lies 

beneath  the  K  shell. 

The  problem  of  collective  oscillations  in  atoms  has 

also  been  considered  by  Wieder  and  Borowitz.     They 

"*"   S.  Wieder  and  S.  Borowitz,  Phys.  Rev.  Letters  l6  (1966)  724, 
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jorioider  the  response  of  an  atom  to  a  l,ong  wavelength 
radiation  field  and  use  first  order  time  dependent  perturba- 
tion theory  to  construct  a  time  dependent  Hartree-Foclc  state 
from  a  stu ._,..- ^y  one.   The  relationship  between  the  indaced 
electron  densi     actuation  and  the  perturbing  radiation 
1    ;  is  linearized  and  a  multipole  expansion  performed. 

The  resonances  in  the  in^i.,--^  iipolu  ...^, .  are  related  to 

the     'ed  otat'j  .  laency  re.     ::es 

are  identified  as  single  Lie  a  high 

frequency  co^..,.^iive  state  is  identified  in  a  r '^l^n  where 
osc'     r  stren  '       ■  single  partic    ■  ■  re 

n   -.   -e.   They  find  t;.     a  col.  -ngths  for 

o   ,  ,  '3      ^   o 


v>^r,.*-^.r^^   Argon  and  Neon  are  207  A,  l46  A  and  16 0  A 

17 
respec     .-      ;.    '  :i   that   if  a   Thomas -Fermi  model 

is  usu       -d  of  :     ;e-Fockj  then  the  collective 

o 
exc~'- — \--^  •"--  —  -  n  is  about  5  A.   In  this  study,  the 

effect  o:     ing  was  not  included  so  we  do  not  know  if  the 

exc      -ns  at  these  frequencies  are  well  defined.   Upon 

comparing  these  values  with  jur^  we  see  that  the  agreement 

is  poor.   Several  reasons  could  account  for  this.   First, 

__ 
'  These  numbers  do  not  appear  in  the  published  work  but 
are  taken  from  S.  Wieder,  Ph.D.  Thesis,  New  York  Univ., 
Oct.  1966.   The  published  paper  reports  a  value  for 
Argon  only.   This  value  is  20  A  which  is  quite  different 
from  the  I96  R   reported  in  the  thesis.   The  discrepancy 
is  attributed  to  the  paper's  being  published  before  the 
completion  of  the  thesis  and  contains  only  rough  estimates. 
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as  we  have  mentioned  our  results  are  not  unique.   It 

could  be  that  had  we  considered  p  state  eigenfunctions 

Instead  of  s  state  ones  we  would  have  found  collective 

modes  in  agreement  with  Wieder  and  Borowitz.   The 

collective  states  of  Wieder  and  Borowitz  are  also  not 

unique  since  they  studied  in  detail  only  the  dipole  terms. 

Higher  multipole  terms  might  have  given  other  additional 

collective  states.   Secondly,  we  did  not  use  the  exact 

Hartree  potential.   Judging  from  the  large  discrepancy 

that  Wieder  and  Borowitz  found  between  the  Hartree-Fock 

solution  and  a  Thomas-Fermi  model  our  approximate  potential 

could  be  a  large  source  of  error.   Finally,  it  could  be 

that  the  collective  state  found  by  Wieder  and  Borowitz 

rapidly  decays  into  single  particle  states.   Our  results 

were  found  by  insisting  that  the  damping  be  small.   If 

we  had  not  demanded  this,  we  could  have  found  solutions 

in  the  same  range  as  Wieder  and  Borowitz. 

There  have  also  been  some  Russian  papers  on  the 

subject.   Unfortionately,  these  papers  are  all  very  sketchy 

and  so  a  detailed  discussion  is  not  possible.   In  one  of 

these  papers,    it  seems  that  the  basic  approach  resembles 

ours  (see  Chapter  V) .   They  use  trial  functions  containing 

36  parameters  and  a  variational  calculation  to  obtain  a 

collective  state  for  Neon  at  68.2  eV.   We  note  that  this 

is  in  very  good  agreement  with  the  first  of  our  values. 
___ 

Y.  Y.  Dmitriev,  L.  N.  Labzowsky,  Physics  Letters  2£A 

(1969)  ^^. 
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It  is  not  clear  from  this  paper  whether,  any  damping  effects 
were  included.   There  is  another  paper  ^  which  alludes  to 
an  experiment  on  inelastic  collisions  of  Argon  ions  in 
which  three  maxima  in  the  cross-section  as  a  function  of 
energy  loss  were  found.   The  author  interprets  these  maxima 
as  evidence  for  collective  states.   Two  of  the  experimental 
maxima  were  40-70  eV  and  220-^20  eV.   The  value  of  the 

third  maximum  is  not  mentioned.   The  second  of  these  is 

o 
in  close  agreement  with  our  first  solution  of  58  A. 


"^^  M.  Ya.  Amusia,  Physics  Letters  l4  (1965)  36 
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APPENDIX   A 


An  Approx ima_t i on  for   the   Density   Matrix 


20 


Suppose  r(x_,p)  represents  a  phase  space  distribution 
function.   Then  the  expectation  value  of  f(x,p)  is 


(A.l) 


<f(x,p)>  ^  r 


j 


f(x,p)  r(x,p)  dx  dp  , 


where  the  hars  signify  operators.   We  can  write 

1  ipq 


(A. 2) 


f(x,p) 


1  ipa 

fq(x)  e      dq  . 


We  note  that 


i  ipq 


e2  ^   fq(x) 


q  S 
e2  ^  f^(x) 


fq(x  +  f)  e 


ipq/2ii 


so  that 

(A. 3)  f(x,p)  =  J  fq(x  +  |)  e  ^   dq 

If  we  take  the  expectation  value  of  f(xjp)  we  have 

i^  ixp 
(A. 4)   <f (x,p)>  =  JJJ  f*{x)    f  (x  +  |)e  ^   e  ^   ^  dx  dp  dq  , 


which  is  the  same  as 


(A. 5)  <f(S,p)>  -   ///  l**(x-  |)  fq(x)  e    e 


^^  i(x-M 


2^ii 


^  dx  dp  dq 


20 


This  approximation  is  based  upon  work  done  by  John  Snygg, 
Ph.D.  ThesiSj  Courant  Inst.  Math.  Sci.^  New  York  Unlv.^1967, 


-152- 


Integrating  over  p  gives 


;a.6) 


<f(x,p)>  ^ 


^  /(x-  §■)  f^(x)  ^{x+   |)  dx  dp 


=  [  f*{x-   |)  f  (x,p)  e    ^   ^(x  +  3.)dx  dp 


J 

We  see  then,  upon  comparing,  (A.l)  with  (A. 6),  that 
the  phase  space  distribution  function  is 

ipq 

ill        (  V  —       -_ 


(A. 7) 


r(x,p)  =  /  e"   ^  fix-   |)  ^(x+  |)  dq  . 


This  expression  is  generally  known  as  the  Wigner  density. 
If    x^  =  X  +  q/2,   Xp  =  X  -  q/2, 

ipCx-^-Xg) 

X-,  +X^  r. 

(A. 8) 


X-,  +Xo 

r(-i^-^  ,p) 


ii 


Tp    (x„)^  (x,  )  dq 


!• 


or 


(A. 9) 


X-,  +x 

r(^^  .p)  = 


iplx-^-x^) 


p(x^,X2)  dq 


where  p(x, ^Xp)  is  the  Dirac  density. 

Hence,  we  can  express  the  Dirac  density  in  terms  of 
the  Wigner  density. 


(A. 10) 


p(x^,X2)  " 


ip(x-, -x„)/ii  x-,+x„ 

^  ^   r(--o-^  .P)  dp  . 


This  expression  is  useful  hecause  it  is  easier  to 
develop  a  closed  form  for  the  Wigner  density  than  for  the 
Dirac  density.   In  order  to  arrive  at  a  closed  form 
expression  for  the  Wigner  density  the  following  identity 
proved  useful*. 
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(A. 11)         r(x,p)  -  I  Tr  (p(x  +  X,  p  +  P)R  ) 

where  the  unbarred  letters  are  c  numbers,  the  barred  letters 
are  operators,   and  R  is  the  reflection  operator. 

We  would  now  like  to  show  that  this  expression  is 
the  same  as  (A. 7).   We  can  rewrite  (A. 11)  as 

-ix-PA   ip-XA         -  r-  P-X  i-  x.p 
(A. 12)   r(x,p)  =  Tr  e        e-       p(x,p)  e   ^     e^     R 

Evaluating  the   trace   gives 

„  -ix-PA   +ip-xA 

(A. 15)       r(x,p)    =   J<x'|e  |x"><x1p(x,p)|x-"> 

-ip-XA   +ix-PA      • 
•<x"'|e  |x^^<x^^1r|x<>   dx'dx"dx"'dxiv 

Now 

-ix-PA   +ip*XA,  -ix'.pA 

<x'le  lx">    =    e  5(x'-x"+X) 

and 

-ip-XA   +ix-PA,  ix'-V-PA 

<x'"|e  lxiv>   =   e  5(x"'-xi^-X) 

and 

<xi^lRlx'>      =      6(x'    +  xiv) 

so    (a. 13)   becomes 

-  —  2x ' - P 
(A.l4)  r(x,p)    =    f  <X+x' 1p|x-x'>    e      ^    ^  dx'     . 


Let   x'    =   2   y« 


-  ryP 


(A. 15)  r(x,p)    =  ^  J  <X-   I  y|p|x+  I  y>    e      ^  dy 

which   is   the   same   as    (A. 7). 
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The  Dirac  density  operator  is  9,(E„-'i)  where  d   is 

2 

the  step  function,  E„  is  the  Fermi  energy,  and  H  ==£■-+  U. 

r  dm. 

Using  an  integral  representation  for  the  step  function  we 
have  -2 

(A. 16)         p(J,p)  =  -  jij-  J  S „.-... dm 

and 

If  U   is    slowly   varying, 

U(x  +X)      =      U(X)    -t  •;)  •  X 

If  we   substitute   this    :  .         [A. 11),   we  have 

(AX8)   1      • t-r  ^.--^'"(vi^ -"('')) 

{IK.LO)      !__..;_        27ri  J      a>fie 


+  Pl£  ^.VU(X)-x) 


1   m  -iiii  m  T^ 

•   ^  Tr  e  R   . 


A   -  7  -x  =  i:  ic^jr  .    ,        B  =  i.^-  +   io.  ^  -f^- 

(A+B)  A      (B+  |-[B,A]+  ^[[B,A],A]) 

(A. 19)  e  -   e      e  '^  ^ 


since   ths  higher  order  conimutators   vanish, 
(A. 20)    [B.A]    =  ^  ^^'-  YjPj    +  lhm2  X:  TjPj 
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(A. 21)  [[B,A],A]  = 


Using  equations  (A.19)j  (A. 20)  and  (A. 21)  to  evaluate  the 
trace  in  (A.l8)  we  have 


(A.  22)  j<  p'le 


[iojS^Y^r^] 


.  -2 
-icbp 

2m 


exp  [VS-  +  i^  21 


in 


+ 


+ 


1  ittD 


..2^.3 


'i^i    2m  ^^  'i  1     6m  ^  'i-^'  ^     -^ 


2   m   '^  'i^i     2m   ^ 


iooS.  Y.  r . 
Since  e    "'"■'"  "^  is  just  a  Taylor  operator  on  the  state 


|p'>,  and  p|-p'>  =  -P'|-P'^^  (A.  22)  is  equal  t 


o 


(A.23)  J  d5p  <p.-.T.lexp  [  '-%-  -  ^^^  -   ^^^-? 


+ 


2  2^2 

CD  Y'Piti  .  Ih  OD-^Y 


2m 


+ 


-I  -P'>  . 


Since  we  have  only  c  numbers  now,  the  only  nonzero 
contribution  is  when 


(A. 2^4) 


p'  -  cuY^  -=  -p'   or  p'  =  -3^  . 


Using  (A. 24) J  the  exponent  in  (A.25)  becomes 

(A. 25) 

V,-    V,    •  r    ih''^a>^lVufx)l^  -. 

which   IS      exp  [ Qir~r^   —  J 


iti'-ar'Yl  i 


2..;5k..w  ni2  _^   > 

since  Y  ^  VU(x) 


If  we  use  correct  5  function  normalization,  (A. 23) 
also  contributes  the  multiplicative  factor,  2^/n-^. 
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Substituting  back  into  (A.18)  gives 


(A. 26)   r(x,p)  =  ^  2^1 


— . —  1  e 
a>i-xe  L 


-i^(E^-  fe  -U(x)) 


F   2m 


iii^otf^  |vu(x)l^/24m 


•  e  ] 

where  we  have  multiplied  by  2  to  account  for  spin. 

Therefore,  for  a  slowly  varying  potential,  the  Dira; 
density  is 


(A. 27)   p(x^,X2)  -  ^  2k 


_d_a^ 
03+1  e 


xp-ix^-Xg) 


i^(E^-  L-  -U(^)+ 


i^jo^ti^  l„,  ,-,|2 


F   2m 


2"^ 


Vu(x) 


where  x 


X^+Xg 
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APPENDIX   B 
The    Evaluation   of   Some    Integrals   which  Appear   In   Method    II 

Consider   the   integrals    over  k,    and   kp,    for   the 


terms    I^,    Ig,    I3  and   l4.  i]  k^-k^l  ^ii^ox  cd^ 

(B.l)      J    d\,     r  d\.  -°' — -P- ^      ^ 

-00  -00 


1   J  "2    /,  2,    2n2       /,  2,_2x2,  2 

(k^+a    )         (kg+P    )    k^ 


CO  00  1 

4Trk^  dk.    J  2(rk^  dk^  J  d   cos    d.^^ 
0  0-1 

I[k^+k2-2k3_.k2]   ^-I^ 


8m(axj_+ri)2)8^^«  p     ^  °°  k-^kg 


-- —  h^ 


dK 


ix^o^^cDg   21  -^Q        ^  ^Q        2    (j^2_^2^2^j^2_^p2^2^2 


1  2  .  ,  2 


"(^-+^278^   (^1+^2)  BETcV^   (1^1 -kg) 

[e  -   8  ] 


T     4-  1       .1  ,1  1  U-+V  ,  u-v 

Let   u   =   k-j^+kg    ,         V  -    ^±-^2    '         1   ^  ~2~  '         2   "^   ~2~ 

The  Jacoblan  =  1/2. 
Also  note  that 
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(B.2)    /  dk3_  J  d.g  -^12 


iCD    CD       ti^ 

2  J    ^^1  J    °^2       ,,  2      2v2/,  2   „2r2"'2 


-00  -00 


(k24<x2)2(k2+p2j\2 


(k^^o-^)    =    (kj_+ia)  (k^-ia). 

We  get 

87r2ni(cD  +03   )    «^  ^  r 
(B.3) -o--— --—  <^du  d\A 


2, 

^i^"2  ^  -OOKD^ 


2      2 
ih  CO.,  cD^u 

J.  id 


4__e__ ^   (u+v)g  p 


(u-v) [u+v+2ia ]2 [u+v-2ia]2 [u-v+2ip]2 [u-v-2ip ]2 

The  integration  over  v  was  performed  by  contour 
integration.   Note  that  there  is  a  singularity  at  v  =  u 
v/hich  corresponds  to  k^  =  0.   The  original  integral  did 
not  have  a  singularity  at  k^  =  Oj  the  singularity  was 
introduced  artificially  hy  the  transformations  made. 
I  therefore  excluded  v  =  u  from  within  the  contour. 
The  result  after  integrating  over  v  is 
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GD 


(B.4)  J  du 


-oo 


2    2 
io)  CO  h   u 

p 

t^Oi^CO^ 

a  p 


(2u-2ia)(4ia)^(2u-2ia+2ip)^(2u-2ia-2i6)^ 


+ 


2ia 


(2u-2ia)^(4ia)^(2u-21a+2ip)^(2u-2ia-2i3)^ 


4ia 


;2u-2ia)(4la)^(2u-2ia+2ip)^(2u-2ia-2ip) 


.^\2 


+ 


4ia 


(2u-2ia)(4la)^(2u-2ia+2ip)^(2u-2ia-2ip)^ 


+ 


4 


xa 


(2u-2ia)  (4ia)  (2u-2ia+2ip)^(2iJ-2ia-2ip) 


.•«^3 


+ 


(-2ip)(2u+2ia+2l6)^(2u+2ip-2ia)^(-4i!3)^ 


+ 


_2u+2ig_ 


(-213)  (2u+2ia+2ip)  (2u+2ip-ia)  (-4ip)' 


4u44ip 


( -2ip ) ( 2u+2ia+2ip )^ ( 2u+2ip-2ia) ^ ( -4ip ) 


4u+4ip 


,2,^..  _...  o,^^3/  ),-•.^2 


(-2i3)(2u+2ia+2ip)^(2u+2ip-2ia)^(-4ip)' 


+ 


_4u44ig^ 


2/  j|^«\3, 


,.  \2 


( -2ip )  ( 2u+2ia+2ip  )  (  -''4ip  ) -^  ( 2u+2i3-2ia) 


-140- 


In  order  to  facilitate  the  integration  over  u^ 
we  decomposed  the  expression  within  the  brackets  into  a 
series  of  partial  fractions.   This  is  done  in  the  follow- 
ing way. 

If  jH — f  i^  3  proper  rational  fraction  (the  degree  of 
the  niAmerator  is  less  than  the  degree  of  the  denominator) 
and  if  f(x)  =  (x-a  )°^(x-b)f^  ...  (x-m)  ,    then 


^f    \  A       A   -,  A-, 

(x-b)P    (x-b)P-l  --^ 


+  ...  + 


M         M  -,  M 

/    ^l^    /    n1^-1  ^  •  •  •    x-m 
(x-m)     (x-m) 


where 


^a-k+l  ~  T^iyi        '      ^(3-k+l  "    (k-l)l  '    "•' 

^V-k+1  ~    Tk^^'l   ' 


and 


1°^       .  .     rt,('^\  /'^._■K^P 


For  example^  consider  the  first  term  in  brackets  in  equation 

(B.-^).  This  is 
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(4ia)^25(u-ia) (u-ia+ip)^(u-ia-ip)' 
^,  (u)    = 


(4ia)^25(u-ia+ip)^(u-ia-ip)' 


V'o(u)    =        - 


(4ia)^25(u-ia) (u-ia-ip) 


2 


^    (u)    = p— ^ p- 

^  (4ia)    2-'(u-ia)(u-ia+ip) 

^  ^  (4ia)      2^p^  2^a   p 


1 


Co   =    ^^(ia+ip)    = ^— ^-^— ^  =   — y-"^---- 

2  ^  (4ia)227iV        2lVp^ 

—  1  "I 

C^   =    ^(ia+ip)    =  -■ ^-^Tf        =   n:0~2~^ 

■^  -^  (4la)    2°^^  2^   a   p 

This  first  term  can  then  he  written  as 

^1        Bg  B^  C^  C^ 

u-ia     (u-ia+ip)      u-ia+ip     (u-ia+ip)      u-ia-i0 

The  integral  over  u  becomes 
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(B.6)       f  du  .--^^—^- e  -L      2 

V 


-00 


h     OO^CUg 


_    ^2  _A_         ^      _-       ^1 

_(u-ia)  u-ia  (u-ia+ip)  u-ia+ip 


B- 


+  — 


+ 


(u-ia+ip)^  (u-ia-ie)  u-ia-i0  (u-'    -'     y 


D- 


+ 


°2  ,__     °1__ 


(u+ia+ip)^         (u+ia+i3)  u+ia+ip 


Tha   contributions   to  ths   coefficients   from  the   various   terms 

are 

term  1 


{iiiaf2^&^ 


-  1 


B«  = 


1 


2        {kxaf2^x^^  2^^aV 


'^V^  ^   2^ OaS"^ 


'2  =    (4i"aVVPp^  ^   2lia\^ 


'l   ^    (4Ta"??"p^  ^   2l0aS'^ 


term  2 


-i 


2   '-   a  29?^ 


,      A^   =    0    ,      B2 


-1 


a  2^"!?^  '      "1 


B.    = 


-5 


a  p52ll    ^ 


^2   = 


-i 


a  2-""p 


m^  ^    ^1  "   ..   ^11«5 


a  2K 
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term  3 


B^  = 


1 


2  "  a22llip5 


^   a  2-^^p 


^?  ~   2  11  ^ 

"^    a  2^  P^ 


^1  ~   2  10  4 


term  4 


-1 


'i  °  ris^ 


B., 


-  1 


^    a  2l0p3 


^2  = 


— '9 — ^ 
a  2^  p 


^1  "  ;;?p5-  "  ^&-p5  =  „2iv   '2 


-X 


C- 


^FS"^  ^  "  a2l2p5 


term  5 


^1  = 


a  2''  p 


5 


^?  "    ^11  T 


%  = 


-5 


a  2l2p5 


'^  "   a   2\^ 


C 


2  ^   a   2^~^ 


c  -  --L 5 

^    a  2^  p5   a   2-^^p 


-11 


12«5  -  «  2l2p5 


term  6 


1 


°2  '  F^^?     °^  °  2^^   p'  a5 


=2'?17J"7     ^I'-^^l-'""' 


2   P"^  a" 


term  7 


-1 


=2  =  -ifpir     %  =  °      =^^?S^ 


B^  =  0 
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term  8 

ot  '      i     r      J,     c 


-1 


^2  -  oiio3^2     ^1  -    ■;^3 


term  9 


term  10 

^2  =  -TTT-   °1  =  °    ^2  ^  ifV-   ^1  =  ° 

Collecting  the  contributions  to  the  coefficients  from 
the  various  terms  we  have 

^1  =  ° 


B,  =  0 


C^  =  0 


^    2^°P^a    "    2-""P"'a 
Note  that 

Cg  =  Dg      C^  -  -  D-^     Ag  =  -  2  Cg- 


*2  = 

a2V 

^2  = 

0 

\ 

=    0 

s  = 

1 

c              ^ 

a  2l°g5 

^2  '   a  210^^ 

°2   = 

i 

^1 

=    0 

ol0o4_ 

and  all  the  rest  are  zero. 
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Carrying  out  the  integrations  indicated  in  equation 

(B.6)^  denoting  —      by  x^,  and  using  the  values  of  the 

1"^  2 
coefficients  just  derived^  we  get  I 

2  2 
-^  ),     f   •  n  1 2  ^      ixti  a 


i>2x 


TRT 


(!-*(«''  /^  )) 


OM  -^2   /  ,Q%     ixh  (a+p) 
2C2iii  X  (a+p)    ^^  ^' 


fm 


(l-*(ii(a4^) 


'-1X 


)) 


2iC,(a-t^)h  X 


52m 


2 


[  -^3^  -7r(a+p)e 


ixh^(a-f^)^ 


m 


(l-*_(ii(a+p) 


"-B5  ))1 


J 


where  *   is  the  error  function! 


*(x)  = 


A 


""   -t^ 

e    dt  , 


0 


¥e  should  notice  hare  the  relatively  simple  form  of 
the  expression,  considering  the  messy  looking  expression 
(B.4)  before  the  integration  over  u  was  done„ 
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APPENDIX  C 

A  Small  Electric  Charge  In  an  Infinlt e  Electron  Oas 

Treated  by  Perturbation  Theory 

ry-     unocc  oc  c 
(G.J)V(r',r%a))  -  f^   [   \_     J_-   u,^(r')  u^Ar")    u^^„(r")  a^„(r') 

.  6   -'i— A_  +  0)) 
—       n 

occ  unocc  E,  ,-E,  II 

+  1;;;  ^u^,(r')a,^,(r")u^„(r")a,^„(r.)  ^^{^^   ^<-)] 

the  u(r)  are  ':      .  nve  functions  for  the 

first  N  levels.   Beyond  that  they  become  continuum  functions 
whsre  we  use  plana  waves  for  the  continuum. 

wher   ;     .3  principal  value. 

Than  the  first  term  of  c,  will  be  of  the  form 

(G.2)g[ii:  Z:;ei^-^  e-^^-^  u^„(r")  a^„(r')  aj-A^-n..) 

4   y^   y~   e^^  ^^  ^         ^^  -^  i  6  (J^— 'L-+  0))] 


Second  term 

^     N_  unocc  ilr'V"    1k"ri     ^ict-Ei," 

(0.5)  ^  [  T:;    F;u,,(r')  a„(r")  e-'^'-  a-^''-  6^{^V-^4<o, 

kT^   k^  +    ti 
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Now  we  want  V*(r',r,a))  r.  J  v(r',r",cD) 
If  V(r',r")  =   7(r'-r"),  then 


r"-r 


dr" 


V(r'-r",cu)  =  r  V(k,a))  e 


^ik<r'-r")  .3 


d^k 


and 


V  (rSr)   =      T(k,cD)  e 


.ik(r'-r")  ^3i,  v(k-)ei^'-(^"-^)d\'d5r' 


V(k,cu)e^^'^'v(k')e"^^''^  6(k-k')  d^k  d\' 


=     J     V(k,a3)  e^^^^'-^)  v(k)  d\  . 

Let  V(r',r",aj)  =  V^  +  V^^  . 

V-p-^  contains  just  plane  waves;  Vj-  contains  the  Coulomb 
functions. 

^ik'(r"  -r) 


*/        .    ...\   _    ['     —       f^,        ^11  „A  A ^3v." 


V  (r^r 


.,0))  .    V   {t',t\o^)   y-^—r   d-^r"  + 
J   X  |r"-r| 


VT.(k,a3)e- 


II 


•v(k)  d^k 


V*   is  the  same  as  the  homogeneous  case. 

We  make  the  simplification  that  there  is  only  one  bound 
state  --  the  ground  state; 

Un  -  2(1-)     e      0 


^0  -  ^^a 


0 


to  evaluate  V*    Then,  the  first  term  of  (C.2)  is 
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,unocG  ,,  .  ,   _,,  .  „  E,,,-E, 


/„  i,\  2f  r-  V- —  ik'.r'   -ik'«r'    /  n  \  -  /  -,  \  .  /  k; '   0   x   1    5  , 


Tha  first  term  of  (C.3)  is 

We  must  integrate 

I r"-r| 
and 

I  r   -r| 


Now 


CO  ^'  -,  r' 


1  '-^^  X/  -,  i 


^-ik'r"cos    e"   ^     ^   (2^+l)(-i)^  j.(k'r")    P   (cos    0") 


^ik"r"    cos    e-   ^     ^      (2i+l)(i)^      J,(k-r")    P,(cos    9") 

i  =  0  '^  ■'^ 

and 


P.     -     Y 


Z     '^      "■  i    /2i+l 
Consider  the   integral 


-1^9- 


(C.6) 


^0  r"-rl 


From  the   angular  part  we   get   a    factor   of     yoT+T  ^  ( -^' --^)'5  ('^-O) 
So    (C.6)   becomes 


(C.7)      f~2    /^     _|^^_(2.,l)(.i)i(|_) 


5/2 


[    r  J,(k'r")   ^.-^  e  °   r"2   dr" 


i' 


^+1 


0 


For  the  time  being  we  will  just  consider  the  i,  =   0 
term  to  see  how  the  problem  develops.   Then  we  have 


,3/2,  z  ,V2.  r  ,.  n...nx  r"°   __^-z^"/ao  _^„2 


2(47r)^/^(f-)^/-[   j.,(k'r") 


0 


0' 


r'"'  dr" 


0 


OD 


,0    -zr"/a, 


+  J  JQ(k'r")  ^--  e  "'  '  °  r"2  dr"] 


Using  jp^(kr)  =  sin  kr  /  kr  and  integrating  by  parts  we  get 


(C.8)   2(47r)5/2^-^) 


5/2 


0 


1  I  il 
ri  k' 


-zr/a„|  a 


;in  k'r-k'cos  k'r 


01.  0 


-zr/a 


,2 

"2 
'O 


ft+k'") 


0 


+ 


— p- [-  —  sin  k'r-k'cos  kr]  + 

^0 
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•zr/a 


0 


+ 


+ 


+ 


[4  +^'']'  ^  '' 

^0 


^z _1 

2(4.)5/2(Z_)5/2 
^0 


[-  T—  COS    k'r  +  k'    sin  k'r] 


-zr/a 


0 


k'[|-  +  k'2j 


z/a„   sin   k'r   +  k'    cos    k'r) 


0 


-zr/a 


0 


Denote  all  this  by  F(kSr)e       .   SO;, 


»         o_  unoGC  .,  ,  ,  -zr/a^   E,  ,-E^ 

(C.9)  v;(r.,r)  =  f  ^  e^'^'-'-'  Uo{r.)  F(k.r)  e     O^jJSl-O^, 

o„  unocc  -zr/a„   -ik"-r'     E„-E,  „ 


5.  h3. 


M.^(co)  :=  J   X(r')  V  (r'.r.co)  cD(r)  d^r  d-'r' 


Try 


-ik^' 
X(r')  =  e      , 


0( r)  =  e 


+ik.r 


Interchanging  the  sum  and  integral, 

ott  ^inocc      „     i(k'-k>r'  ^ 

(CIO)    H^    (a.)   =  g-     ^   (   J    e  UQ(r«)    d^'>) 


J  F(k'r) 


-zr/aQ     ik-r  ,^1-^0 


-i(k+k">r»      , 
k 
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-zr/a^  ik-r         ^n"\" 
F(k"r)  e     ^  e    d^r]  ^^[-\^-   +^1 


Consider  the  integral  over  r'. 


®0 


3/2 


-zr'/an   "  l(k'-k)T' 
e      ^  e  ^^^' 


d^r' 


i(k'-k)'r' 


u„(r')  d-'^r  . 


Ths  integral  when  evaluated  is 


(C.ll)        47r(^)5/2 


^°'     [-2  -(k'-k)2+  4(k'-k)2] 
^0 


and 


(C.12)    u„(r')  e 


-i(k+k").r"   . 


0' 


d^r' 


^•^(^)5/' 


4 


'0 


[£   -(k+k")2]2+4(k+k")2 


'0 


To  estimate  the  order  of  magnitude  of  this  term  we  note  that 


k'. 
mm 


k„  ~  10 

F 


00 


0 


and 


k'  »  k  . 


Therefore^  this  term  is  of  the  order 

,V2 


^ ^ r^ 


(k-)- 


?r 


k' 


-5/2 


and  the  term  is  a  maximum  for  k'  =  k„.  So,  the  maximum  value 
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-S/2 
is  about  kp-^^  .   Consider 


-zr/a^   ik-r   -^ 
F(k«r)  e       e     d-^r 


a^  ~  10     and   F(k'r)  is  sinusoidal. 

Therefore  the  integrand  must  be  very  small  unless 

-8 
r  is  small^   r  <  10 

Also,  judging  from  results  for  the  homogeneous  gas, 

it  seems  reasonable  to  assume  k  small,  i.e.  k/k^  <<  1, 

Therefore 

ik*r 
e      -^^  1  . 


This  would  justify  considering  only  the  i,   -^  0  case. 

o 

The  minimum  value  of  k'  is  k^^^^^^  ~  10  .   Also  z/Bq  ~  1' 
¥e  note  that  the  first  and  last  terms  in  F(k'r)  cancel. 


(0.15)      F(k-r)    =   2(47r)5/2(^)5/2 


1  1 


^0 

^2 
r    -  — ^—  sin   k'r   -   —  cos    k'r   -  -^os   k'r+k'sin   k'r  +  ~] 


2(4.)5/2(|-)V2 

--^- ^ \   [sin  k-r[k'    -  -^]-  f-  cos   k'r+  f -] . 


0    pQD  -zr/aQ  2 

Evaluating  F(k'r)    e  47r  r     dr,    we  have: 

^0 
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(0.14) 


2(4.)3/2(^.)3/2 


tL 


0 


2 

;o 

2 


-   8v  |-[  ^2-  -k'^] 


/   [   -o-  +k'"] 


'0 


.2i2 


0 


a;;k' 


Stt 


+  z/a 


0 


J 


oo 


-zr/a 


0 


F(k"r)    e  ^  -^lir  r     dr      is    the   same   with   k'    replaced   by   k". 

^  Now  we   have   to   sura  over  k"    and   k"   with  both  k'    and   k" 

o 

going  from  k„   to   oo ,    and   k^  -%.  10        for  a   high  density   gas. 
Wewill   therefore   approximate    (C.l4)    in  order  to  get   an 
estimate   of   the  magnitude. 

The   expression   in   the   brackets    is   ~  lA'    ^or  k'    very 
large.       (The  maximum   of   the   expression   is   for  k'    =   k^  ^  z/Qq. 
Therefore 


00 

J  F(k'r) 
0 


-zr/a 


0 


hir   r     dr  ~  — 


^0 
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k'"   k' 


k.5/2  -7/2 

Tr5 "  =  ^F 


=     maximum  value. 
The  magnitude   of   the  perturbation   terin   is 

,-5/2  .  ^,-7/2,  f-^.-6 


00 

k'=kp 


and  the  maximum  value  is  (10   )   =10 

Therefore  for  z  positive  there  is  no  measurable  effect 
The  inclusion  of  additional  terms  would  only  give  more 
terms  of  this  order  of  magnitude. 
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